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SUMMARY 


A Navy  T-2  jet  trainer  aircraft  vaa  instrumented  to  measure  and  record  all 
airframe  motion  variable?.  Motion  time  histories  were  recorded  for  small 
perturbation  responses  to  a variety  of  carefully  selected  pilot  inputs.  Two 
flight  conditions  were  considered:  a relatively  high  speed  power  approach 

configuration  and  a high  speed  cruise  condition.  Noise  levels  in  the  data  were 
found  to  be  significant,  and  it  was  determined  that  measurement  noise  was  a 
larger  problem  in  this  case  than  process  noise.  A unique  noise  problem  was 
the  high  level  of  measurement  noise  on  the  control  Inputs. 


Although  all  motion  variables  were  recorded,  only  longitudinal  swtlon  time 
histories  were  used  for  analysis.  Three  digital  computer  parameter  Identifi- 
cation techniques  were  applied  to  the  data:  modified  Newton-Raphson  (in-house 
at  NAVA IRDEVCEN ) , Kalman  filtering/smoothing  (Calspan  Corp.) , and  maximum 
likelihood  (Systems Control,  Inc.).  The  objective  was  the  determination  of  the 
relative  accuracy  and  utility  of  the  three  techniques.  Volume  I reports  the 
data  gathering  effort  and  the  modified  Newton-Raphson  analysis. 


For  the  in-house  analysis,  a linear  ssathematical  mode]  was  used, 
containing  nine  stability  derivatives  to  be  identified.  In  general,  the 
identification  effort  was  successful  in  matching  model  response  to  flight 
data  and  obtaining  consistent  sets  of  stability  derivatives  from  the  various 
data  segments,  but  the  time  history  matches  are  far  from  perfect  and  some 
anomaliej  do  exist  in  the  stability  derivative  results. 
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INTRODUCTION 

The  problem  of  extracting  aerodynamic  parameters,  such  as  stability 
derivatives,  from  aircraft  flight  data  has  been  the  subject  of  considerable 
study  over  the  years.  In  the  past  decade,  however,  the  use  of  digital  computer 
parameter  identification  techniques  has  given  new  impetus  to  the  technology. 

The  state  of  the  art,  except  for  the  most  recent  developments,  is  summarized 
in  reference  (a) . 

The  primary  objective  of  this  program  was  the  assessment  of  the  relative 
accuracy  and  utility  of  several  existing  parameter  identification  techniques. 
Secondary  objectives  were  the  establishment  of  one  consensus  set  of  stability 
derivatives  for  the  T-2  aircraft,  development  of  one  identification  technique 
for  continuing  in-house  use  at  NAVAIRDEVCEN,  and  the  establishment  of  a 
carefully  and  completely  documented  file  of  flight  data  to  be  made  available  to 
anyone  pursuing  parameter  Identification  technique  development.  Although  not 
part  of  the  original  objectives,  the  mri«t  significant  result  of  this  current 
work  may  be  the  development  of  NAVAIRDEVuEN  expertise  (in  data  acquisition  and 
use  of  identification  techniques)  which  will  be  essential  in  the  current 
effort  on  development  of  high  angle  of  attack  (stall/post-atall/spin) 
parameter  identification  processes. 

It  should  be  emphasized  that  no  attempt  was  planned  or  made  to  originate 
new  identification  techniques  or  even  substantially  modify  existing  ones. 

Rather,  existing  techniques  were  used,  modified  only  slightly  for  adaptation 
to  this  particular  application.  In  recent  years,  the  opinion  that  the  low 
angle  of  attack  aircraft  identification  problem  had  been  "solved"  has  been 
expressed  by  numerous  investigators  in  various  forusis.  Thus  there  seemed  to 
be  no  need  to  attempt  to  develop  any  new  identification  techniques  for  this 
program. 

Since  a portion  of  the  program  objective  was  the  establishment  of  an 
in-house  identification  capability  for  the  future,  the  existing  telemetry  and 
data  recording  equipment  at  NAVAIRDEVCEN  was  used  for  the  project,  although 
this  equipment  was  by  no  means  state-of-the-art.  This  constraint  conflicted 
with  the  desire  of  producing  carefully  controlled,  low  noise  level  flight 
data  • a type  of  flight  data  which  has  been  essentially  unavailable  to  many 
investigators  in  the  parameter  identification  field.  The  most  significant 
problem  with  generally  available  flight  data  in  the  past,  however,  has  been 
a lack  of  documentation.  For  this  project  all  relevant  information  pertaining 
to  flight  circumstances  and  instrumentation  has  been  recorded. 

The  need  for  an  in-house  capability  in  parameter  identification  perhaps 
merits  further  comment.  Although  it  has  not  been  active  in  in-house  experimental 
flight  dynamics  research  programs  prior  to  the  effort  reported  here,  NAVAIR- 
DEVCEN is  designated  (reference  (b))  as  the  lead  laboratory  for  flight  dynamics 
research  and  development  work  in  the  Navy.  Consequently,  it  would  be  appropriate 
and  desirable  for  NAVAIRDEVCEN  to  have  a capability  for  in-house  experimental 
flight  dynamics  work,  and  parameter  identification  would  probably  be  required 
in  such  work. 
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In  particular,  NAVAIRDEVCEN  has  planned  to  construct  a variable  stability 
research  aircraft  and  prototype  trainer  based  on  the  T-2.  A parameter 
identification  capability  is  necessary  for  the  calibration  and  verification 
of  such  a simulator.  A consensus  set  of  T-2  stability  derivatives  as  a 
function  of  flight  condition  would  be  necessary  for  the  programming  of  the 
variable  stability  control  system,  so  two  portions  of  the  results  of  this 
work  have  application  to  the  variable  stability  project. 

J 

Volume  I of  this  report  concerns  the  acquisition  of  the  flight  data  and 
the  in-house  identification  effort,  while  Volume  II  reports  on  the  contractor 
identification  efforts. 


SELECTION  OF  IDENTIFICATION 
TECHNIQUES 

It  is  generally  recognized  that  modern,  digital  computer.  Identification 
techniques  can  be  divided  into  two  types:  maximum  likelihood  and  Kalman 

filtering.  Reference  (c)  provided  an  explanation  of  this  division.  Within 
these  two  broad  divisions  there  exist  many  techniques  which  are  somewhat 
different  in  specifics  although  similar  in  principle,  and  in  fact  there  are 
numerous  techniques  which  are  quite  similar  although  referred  to  by  different 
names.  All  these  techniques  are  "output  error"  types;  equation  of  motion 
error  techniques  and  non-digital  (l.e«,  analog  matching)  techniques  were  not 
considered  in  this  study. 

Probably  the  most  popular  of  the  digital  techniques  today  is  the  modified 
Newton-Raphson  technique,  which  is  a member  of  the  maximum  likelihood  "family", 
and  in  fact  becomes  a maximum  likelihood  method  if  a properly  calculated 
weighting  matrix  is  used.  Considering  the  relatively  straightforward 
nature  and  the  wide  acceptance  of  the  modified  Newton-Raphson  technique,  it 
was  decided  to  use  this  technique  in-house  at  NAVAIRDEVCEN. 

Processing  of  the  data  by  the  most  advanced  identification  technique  at 
NASA -Ames  Research  Center,  NASA -Langley  Research  Center,  and  NASA -Flight 
Research  Center  was  considered  highly  desirable,  but  pressures  of  other  work 
and  limited  manpower,  and  delays  in  the  gathering  of  the  data  by  NAVAIRDEVCEN, 
made  NASA  participation  uncertain.  To  ensure  that  at  least  three  different 
identification  techniques  would  be  used,  two  contractors  were  selected  to 
perform  analysis  concurrently  with  the  in-house  effort. 

Systems  Control,  Inc.,  of  Palo  Alto,  California,  was  contracted  to  use 
their  version  of  maximum  likelihood  (reference  (d)>;  Calspan  Corp.  of  Buffalo, 
New  York,  was  paid  to  apply  their  Kalman  filtering/smoothing  (reference  (e)) 
method  to  the  data. 

It  was  decided  that  the  in-house  analysis  effort  would  have  a different 
philosophy  than  the  contracted  efforts.  In-house,  the  data  would  be  run 
through  the  identification  procedure  without  any  substantial  operator  inter- 


face - the  program  would  simply  be  "on  its  own".  The  contractors,  however. 
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and  to  use  any  level  of  human  involvement  which  they  felt  appropriate.  Major 
modifications  to  the  programs  were  rendered  impossible  by  cost  and  time 
constraints,  in  addition  to  being  contrary  to  the  philosophy  of  the  program. 


EXPERIMENTAL  EQUIPMENT  AND 
CONDITIONS 


THE  AIRCRAFT 

The  flight  data  used  in  this  report  were  obtained  with  YT-2B  BuNo  144218, 
shown  in  Figures  1 and  2.  The  aircraft  was  constructed  by  Rockwell  International, 
Columbus  Division  (formerly  North  American) , as  a prototype  of  the  model  T-2B, 
and  is  essentially  identical  to  the  T-2B  and  T-2C  with  respect  to  aerodynaad.es. 

The  geometric  parameters  of  the  aircraft  are  given  in  Table  I. 

The  most  significant  features  of  the  aircraft  configuration  are  the  straight, 
mid-fuselage  wing  and  the  conventional  elevator  mounted  on  a fixed  cruciform 
tall.  The  flaps  are  send -Fowler  single-slotted. 

Naval  Air  Facility,  Warminster,  was  the  base  of  the  aircraft  throughout 
the  flight  program,  and  the  majority  of  flights  were  conducted  in  the  south- 
eastern Pennsylvania  and  southern  New  Jersey  areas. 

INSTRUMENTATION 

A complete  set  of  motion  transducers  was  installed  in  the  YT-2B  for  the 
sole  purpose  of  obtaining  flight  time  histories  for  parameter  identification 
analysir.  The  locations  and  serial  numbers  of  the  transducers,  as  originally 
installed,  are  given  in  Table  II.  The  measurement  ranges  of  the  transducers 
are  presented  in  Table  III.  As  of  10  May  1973,  the  changes  listed  in  Table  IV 
became  effective. 

In  addition  to  the  transducers  listed  in  the  tables,  the  aircraft  was 
equipped  for  the  measurement  of  airspeed,  altitude,  angle  of  attack,  sideslip, 
and  control  positions.  For  the  purposes  of  this  project,  the  YT-2B  was  re- 
equipped with  the  noseboom  constructed  by  North  Asmrlcan  at  the  time  of  the 
original  flight  testing.  The  boom  included  a pitot-static  system  and  * and 
P vanes. 

The  pressure  altitude  transducer  was  Intended  to  measure  perturbation 
altitudes  ±1200  feet  about  a reference  altitude,  which  was  "locked  in"  when 
the  pilot  selected  the  "altitude  reference  engage"  control.  The  transducer 
malfunctioned,  however,  and  would  measure  only  +1200,  -600  feet  from  reference. 

A reference-and-cllmb  procedure  was  used  to  make  the  effective  range  ±900  feet. 

Airspeed  was  also  calculated  from  the  dynamic  pressure  transducer  at  the  tip 
of  the  noseboom.  Two  independent  transducers  were  used  for  "coarse"  and  "fine" 
measurements.  Coarse  airspeed  had  a range  of  0-500  kt,  while  0-50  kt  was 
the  range  for  fine  airspeed  which  was  automatically  reset  to  zero  at  each 
50  kt  upper  limit,  thus  describing  a "sawtooth"  function.  This  resetting 
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'ICURE  1.  YT-2P.  V.uNo  144218,  CLEAN  CONFIGURATION 


FIGURE 
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TABLE  I 

GEOHETRIC  PARAMETERS  OF  FULL-SCALE  T-2 
(from  reference  (f;) 


WING 

SW 

bw 

ARV 
\ w 
rw 


c 


r 


c 


t 


c 


w 


i 


w 


Total  area  (includes  flap,  aileron  and 
39.39  ft2  covered  by  fuselage) 

Net  surface  area  (wetted) 

Span  (perpendicular  to  plane  of 
symmetry)  including  tiptanks 

Aspect  Ratio 

Taper  Ratio 
Dihedral  Angle 

Chord  (in  streamline  direction) 

Root  (Wing  Sta.  0) 

Tip  Chord  (Wing  Sta.  214.242) 

(Equivalent) 

Mean  aerodynamic  chord 

(Wing  Sta.  95.078) 

Location  of  257.  MAC 
Sweepback  of  257.  element 

Incidence  angle 

Root  Chord 
Tip  Chord 

Airfoil  Section  (root  and  tip  in 
streamline  direction) 

*NAA  Modified 


254.86  ft2 
424.85  ft2 

38.13  ft 

5.07 

.496 

+3* 

114.20  in 
56.63  in 

88.88  in 

F.S.  219.697 
2*17' 


2° 


-1° 

NASA64]A212 
a - .8* (MOD) 
(flaps  and 
ailerons  rigged 
3°  up) 
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TABLE  I (CON'T) 

GEOMETRIC  PARAMETERS  OF  FULL-SCALE  T-2 
(from  reference  (£)) 


Rate  of  Taper 

FIAP  (Data  for  One) 

Type 

Area 

Span  (perpendicular  to  plane  of  syometry) 

c^  Inboard  chord  (Wing  Sta.  27.09) 

c Outboard  chord  (Wing  Sta.  127.54) 

o 

Cj/cw  Ratio  flap  chord  to  wing  chord  (avg.) 
b^/bw  Ratio  flap  span  to  wing  semi-span 

6^  Flap  deflection,  maximum  (from  uprigged 

position) 

Flap  in  neutral  position 

AILERON 


Type 

S Area  (aft  of  hinge  line  and  including  tab) 

s 

b Span  (perpendicular  to  plane  of  syninetry) 

c^  Inboard  chord  (Wing  Sta.  T.28.69) 

cq  Outboard  chord  (Wing  Sta.  208.26) 

ca/cw  Ratio  aileron  chord  (aft  H.L.)  to  wing  chord 

ba/bw  Ratio  aileron  span  to  wing  semi -span 
2 


0.2671 

Single  Slotted 
22.78  ft2 

101.75  in 

39.39  in 

29.63  in 

.37 

.475 

33° 

3®  Up 

Straight  Sided 
9.5  ft2 

79.57  in 

20  in 

14.66  in 

.25 

.374 
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' Aileron  deflection,  maximum  (from  neutral 

position) 

Aileron  in  neutral  position 
Aerodynamic  Balance 

Sb  Balance  area  forward  of  the  H.L. 

(including  507.  of  fabric  seal) 

c, /c  Ratio  balance  chord  to  aileron  chord 

b a 

Static  balance 

Irreversible  full  power  system 
AILERON  TRIM  TAB 

Ground  adjustable  fixed  tab  on  each 
aileron 

S Area  (each) 

31 


-12*  Up,  +13*  Dn 
3*  Up 

Sealed  paddle  balance 

4.45  ft2 

.42 

Weighted  paddle  balance 
Hydraulic 

.07  ft2 


HORIZONTAL  TAIL* 


sh 

Total  area  (Includes  3.07  ft2  covered 
by  vertical  tail  and  fairing) 

72.29  ft2 

sneth 

Net  area 

146.38  ft2 

Ah 

Net  surface  area  (wetted) 

146.38  ft2 

bh 

Span 

17.91  ft 

ARh 

Aspect  Ratio 

4.42 

Xh 

Taper  Ratio 

0.50 

’ h 

Dihedral  angle 

0* 

Ah 

Sweepback  of  257.  clement 

15* 

♦Percent  lines  base  on  horizontal  prior  to  addition  of  trailing  edge  extension. 
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TABLE  I (CON'T) 

GEOMETRIC  PARAMETERS  OF  FULL-SCALE  T-2 
(llrom  reference  (£)) 


Chord  (in  streamline  direction) 

Root  (H.T.  Sta.  0) 

Equivalent  tip  chord 

(H.T.  Sta.  106.483) 

Mean  aerodynamic  chord 

(H.T.  Sta.  47.78) 

Incidence  angle 

Airfoil  section  (root  and  tip  in  streamline 
direction) 

Tail  length  (.25  c to  .25  c.  ) 

w n 

HORIZONTAL  STABILIZER 

S Area  stabilizer,  total 

s 

i Stabilizer  incidence  angle 

8 

ELEVATOR 

S Total  area  (excluding  balance  area  forward 

of  the  hinge  line) 

b Span  (between  equivalent  chords) 

(one  elevator  only) 

c^  Inboard  chord  (B.P.  3.906) 

cq  Outboard  chord  (B.P.  105.877) 

cc/ch  Ratio  elevator  chord  (aft  H.L.)  to 

horizontal  tail  chord 

b /b.  Ratio  elevator  span  to  horizontal  tail 

c n 

span 


ch 

*h 


64.61  in 
33.05  in 

50.447  in 

0* 

NASA  65A012 
202.58  in 


42.5  ft2 
0° 

21.00  ft2 

101.97  in 
18.85  in 
10.52  in 

.310 

0.936 
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TABLE  I (CON'T) 

GEOMETRIC  PARAMETERS  OF  FULL-SCALE  T-2 

fie 

(from  reference  (f)) 
Elevator  deflection  maximum 

27*  Up,  15*  Dn 

Boost:  Push  force  2.95:1 

Hydraulic 

Pull  force  2.95:1  to  8 lbs 
then  6.0:1 
Static  balance 

Weighted  Leading  Edge 

Aerodynamic  balance 

1 i 

Overhang 

Sb 

Balance  area  forward  of  hinge  line 

5.72  ft2 

\ j 

cb/ce 

Ratio  balance  chord  to  elevator  chord 

0.322 

Nose  factor 

i 

0.60 

Point  of  tangency  for  nose  factor  is 
at  elevator  hinge  line 

ELEVATOR  TRIM  TAB 

St 

Area  (each) 

2.36  ft2 

bt 

Span,  Equivalent  (B.P,  8.93  to  54.53) 

46.10  in 

1 ] 

Ct 

Chord,  constant 

6.5  in 

\/be 

Ratio  tab  span  to  elevator  span 

0.462  in 

K 

Tab  deflection 

L.H.  10°  Up, 

VERTICAL  TAIL 

13°  Dn 
R.H.  0°  Up, 
13*  Dn 

S 

Total  area  (includes  4.38  ft2  blanketed 

40.33  ft2 

V 

by  fuselage  plus  2.14  ft2  blanketed  by 

^nctv 

horizontal  tail) 
Net  area 

33.86  ft2 
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TABLE  I (CON'T) 

flECMETRIC  PARAMETERS  CF  FULL-SCALE  T-2 
“ (from  reference  ( OJ 


Net  surface  area  (wetted) 

Net  surface  area  of  dorsal  fin  (wetted) 
Span,  unblanketed 
Aspect  ratio 
Taper  ratio 

Chord  (in  streamline  direction) 

Root  (W.P.  + 33.000) 

Equivalent  Tip  Chord  (W.P.  +129.41) 

Mean  aerodynamic  chord  (W.P.  + 73.92) 

Sweepback  (257.  chord) 

Airfoil  Section 

Tail  length  (.25  cw  to  .25  cy) 


VERTICAL  FIN 
S* 


RUDDER 

S 


Area  (including  2.14  ft2  blanketed  by 
horizontal  tail  and  excluding  dorsal  fin) 

Angle  with  respect  to  airplane  plane  of 
symmetry 


Total  area 


STu  Upper  surface 


79.18  ft2 
18.12  ft2 
8.04  ft 
1.80 
.375 

78.14  in 
29.38  in 

58.47  in 

30° 

NASA  63A012 
194.05  in 

29.87 

0* 

9.13  ft2 
3.23  ft2 
5.90  ft2 


Sri  Lower  surface 
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TABIZ  I (CON'T) 

GEOMETRIC  PARAMETERS  OF  FULL-SCALE  T-2 
~(ir<m  reference  (£)) 


Span,  equivalent 


bru  Upper  surface 
bri  T,ower  surface 

cru  Upper  chord  (W*P.  96.00) 

cri  Lover  chord  <W.P.  + 9.91) 


c /c  Ratio  rudder  chord  (aft  H.L.)  to  vertical 

r v tail  chord 

cf/cv  Upper  surface  @ W.P.  96.00 

c /c  Lower  surface 
r'  v 

Rudder  deflection,  maximum 
Boost 

Aerodynamic  balance 

Balance  area  forward  of  hinge  line 

c./c  Ratio  balance  chord  to  rudder  chord 

b r 

cb/cru  Upper  surface  @ W.P.  96.00 


c./c  Lower  surface 
t>  r 

Static  balance 
Nose  factor 

Point  of  tangency  for  nose  factor  is  at 
rudder  hinge  line 

RUDDER  TRIM  TAB 
St  Area 


31. 94  in 
42.99  in 
12.59  in 
22.45  in 

.266 

.250 

25*  Rt,  25*  Lt 

None 
Overhang 
2.41  ft2 

.234 

.24 

Weighted  leading  edge 
0.40 


1.60  ft2 


1 

1 
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TABLE  I (CON*T) 

GEOMETRIC  PARAMETERS  OF  FULL-SCALE  T-2 
(from  reference  (£)) 


b/b 
t r 


FUSELAGE 


w. 


L/D 

CANOPY 


r 

C 


A 

c 


NACELLE 


Span,  equivalent  (W.P.  14.94  to  W.P.  53.00) 

Chord,  constant 

Ratio  tab  span  to  rudder  span 

Tab  deflection,  maximum 

Length  (actual) 

Maximum  frontal  area  (basic  fuselage) 
Maximum  width  (basic  fuselage)  F.S.  169 
Maximum  depth 

Basic  fuselage  over  canopy  (F.S.  169) 
Including  ducts  (F.S.  214) 

Net  surface  area 

Fineness  retio  (actual) 

Length  (actual) 

Maximum  frontal  area 
Net  surface  area 
Fineness  ratio  (actual) 


1 

n 

V 
* V 


Length  (actual) 
Maximum  frontal  area 


38.06  in 
6.0  in 
.508 

7°  Rt,  7“  Lt 

34.58  ft 
15.75  ft2 
54  in 


88.1  in 

73.9  in 
221.11  ft2 

5.91 

19.75  ft 
3.70  ft2 

73.10  ft2 

8.8 

23.71  ft 
10.50  ft2 
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TABLE  I (CON*T) 

GEOMETRIC  PARAMETERS  OF  FULL-SCALE  T-2 
(from  reference  (f)) 


A 

n 

Net  surface  area 

206.0  ft2 

Inlet  area  (includes  gutters) 

3.1  ft2 

L/D 

n 

Fineness  ratio  (actual) 

5.025 

SPEED  BRAKE  (Data  for  one  side  only) 

Type 

One  Piece 

Location 

Side  of  Aft  Fuselage 

Number 

Two 

sj 

Area  (Planform) 

8.00  ft2 

pj 

Area  (frontal) 

4.24  ft2 

Maximum  deflection 

32°* 

TIP  TANK 

(Data  for  one  tank  only) 

ht 

Overall  length 

142.75  in 

dtt 

Maximum  diameter  (Tank  Sta.  61.875) 

20.00  in 

L/D 

Fineness  ratio 

7.14 

Sstt 

Side  area  (projected) 

14.1  ft2 

SPtt 

Planform  area  (projected) 

14.2  ft2 

Volume 

15.3  ft2 

Att 

Total  Surface  Area 

44.30  ft2 

^nettt 

Net  Surface  Area  (wetted) 

42.40  ft2 

*Noto:  Manufacturers  specification  is  32°.  By  actual  measurement,  test  A/C 

maximum  speedbrake  deflection  is  50*.  Reference  (g)  states  that  AFC  103 
ordered  extension  of  speed  brake  deflection  to  50°  for  production  T-2  aircraft 
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TABLE  II 

INSTRUMENT  LOCATIONS  ON  YT-2B  FOR  AIRFRAME  DYNAMICS  IDENTIFICATION  PROGRAM 

(prior  to  10  May  73) 


Function 

Name 

Fuselage 

Station 

Wing 

Station 

Water 

Line 

n 

X 

Greenleaf  Serial 
1655000-A  #674 

37. 

-4. 

-11. 

nz 

(c.g.) 

Giannini  #1816 
24117P-3.5-20 

207. 

-1. 

• 

H 

Csl 

nz 

• 

nose  (for  q) 

Giannini  #1819 
24117P-3.5-20 

34. 

-6. 

-11. 

nz 

• 

tail  (for  q) 

Giannini  #1837 
24117P-3.5-20 

366.5 

0 

+9. 

• 

q 

Stahara  Ang.  Acc. 
AA196-8-350 

230. 

0 

-23.5 

8 

Giannini  Model 
Serial  5812  3416 

56. 

0 

5. 

ti 

Giannini  Model 
Serial  3812  3416 

56. 

0 

5. 

r 

Gyro  Dynamics 
Serial  #62 

47. 

0 

3.5 

P 

Norden  RG  228 
Serial  #185 

52. 

0 

2.5 

q 

Norden  RG  228 
Serial  #183 

196. 

0 

-22.5 

ny 

• 

nose  (for  r) 

Edelcliff  Model  7-30 
Serial  #367 

31.5 

-6. 

-11. 

n 

y 

tail  (for  r) 

Edelcliff  Model  7-30 
Serial  #366.5 

366.5 

0 

+9. 

n 

y 

(c.g.) 

Edelcliff  Model  7-30 
Serial  #4444 

204. 

-1.5 

-19. 

nz 

• 

wing  tip  (for  p) 

Giannini  #1685 
24117P-3.5-20 

249. 

206.5 

15. 

nz 

• 

wing  tip  (for  p) 

Giannini  #1860 
24117P-3.5-20 

249. 

-206.5 

15. 

NOTE;  Negative  wing  stations  are  toward  left  wing.  Negative  water  line 
positions  are  below  fuselage  reference  line. 
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TABLE  III 

MEASUREMENT  RANGES  FOR  YT-2B  TRANSDUCERS 


Measurement 


Range 


n 


z 


±3.5g 


q 

8 


r 

q 

p 


±1.0g 
±1.0g 

±20  deg/so'* 

±15  deg 
±20  deg/flec 
±450  deg/aec2 
±90  deg 

±20  deg/sec  (prior  to  5/10/73) 
±45  deg/sec  (effective  5/10/73) 


25 


NA  DC-74181-30 


function  developed  problem , however,  end  ceueed  a one  noise  in  the  velocity 
signal  at  the  switching  points.  A position  error  correction  for  the  airspeed 
measureaent  was  supplied  by  the  aanufaeturer  in  reference  (h). 

The  angle  of  attack  and  sideslip  vanes  were  located  23.1  feet  and  23.4  feet, 
respectively,  forward  of  the  noalnal  center-of -gravity.  The  ranges  of  the 
measurements  from  these  vanes  were  adjustable,  but  were  set  at  ±10*  for  both 
vanes  throughout  this  flight  program. 

While  all  gyros  and  accelerometers  were  aligned  with  respect  to  the 
fuselage  reference  line  indicated  by  the  aircraft  manufacturer,  the  n vane 
was  referenced  to  the  nosebooa.  It  was  calculated  that  the  noseboom  Itself 
was  inclined  -4.0*  relative  to  the  wing  chord  line,  from  which  true  angle  of 
attack  should  be  measured . 

Aileron,  rudder,  elevator,  flap,  and  speedbrake  positions  were  measured 
via  potentiometers  covering  the  full  range  of  available  movement. 

Note  that  pitch  angular  acceleration  was  measured  twice:  once  directly 

with  the  Statham  angular  accelerometer  and  once  via  differential  normal 
acceleration  between  nose-mounted  and  tail-mounted  normal  accelerometers. 

TELEMETRY,  RECORDING,  AND  DATA  REDUCTION 

Transducer  signals  were  simultaneously  recorded  onboard  the  aircraft  and 
broadcast  to  the  ground  station.  The  telemetry  system  was  of  the  FM-FM  type 
and  the  onboard  recorder  was  an  Astro-Science  Corporation  MARS  2000.  Since 
only  12  discrete  IRIG-standard  subcarrier  frequencies  were  available,  many 
of  the  channels  had  to  be  assigned  to  a pulse  amplitude  modulation  system  and 
then  sampled  at  30  samples/sec  by  a commutator,  with  the  coanaitator  output 
modulated  on  the  70  KHZ  subcarrier.  Table  V lists  channel  assignments.  The 
various  subcarrier  signals  were  then  modulated  together  and  recorded  as  a 
single  channel  on  the  tape. 

Demodulation  and  analog-tc-dlgital  conversion  could  not  be  accooplished 
at  NAVAIRDEVCEN  without  procurement  of  additional  equipment,  so  these  operations 
were  performed  at  Naval  Air  Test  Center.  The  digital  data  tapes  were  then 
provided  as  input  to  the  NAVAIRDEVCEN  digital  computer  facility. 

A digital  cosqtuter  program  for  data  reduction  purposes  was  developed. 

This  program  transformed  the  digital  information  into  CDC  6600  format,  corrected 
for  aircraft  instrumentation  offsets,  decallbrated  the  data  into  physical 
units,  and  referenced  all  motion  variables  to  trim  values.  A printout  of  this 
program,  configured  for  a typical  run,  is  given  as  Appendix  A. 

The  calibration  constants  used  in  this  process  for  the  gyros  and  accelerometers 
were  developed  from  tests  at  NAVAIRDEVCEN  using  rate  tables  and  similar 
equipment.  Every  transducer  was  tested  independently.  Control  surfaces  and 
c and  P vanes  were  calibrated  using  an  inclinometer.  Airspeed  and  altitude 
transducers  were  calibrated  in  accordance  with  procedures  recommended  by 
manufacturers  using  pressurized  air  sources. 
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TABLE  V 

CHANNEL  ASSIGNMENTS  FOR  YT-2B  TELEMETRY  AND  RECORDING 


Channel 


FM  10 
FM  14 
FM  22 
FM  40 
FM  70 


.3KHZ 

.7KHZ 

.3KHZ 

.OKHZ 

. 9KHZ 

.4K3IZ 

. 35KHZ 

.5KHZ 

.SKHZ 

.OKHZ 

.OKHZ 

.OKHZ 


Assigned  Variable 

angle  of  attack 
sideslip  angle 
heading 
pitch  angle 
pitch  rate 
roll  rate 
roll  angle 

normal  acceleration  (port  wingtip) 

pitch  angular  acceleration 

normal  acceleration  (starboard  wingtip) 

aileron  position 

PAM 


PAM1 

PAM2 

PAM3 

PAM4 

PAM5 

PAM6 

PAM7 

PAM8 

PAM9 

PAM10 

PAM11 

PAM12 

PAM13 

PAM14 

PAM15 

PAM16 

PAM17 

PAM18 


zero  reference 
full  scale  reference 
flap  position 
speed  brake  position 
flight  path  acceleration 
side  acceleration  (eg) 
normal  acceleration  (eg) 
normal  acceleration  (nose) 
normal  acceleration  (tail) 
side  acceleration  (nose) 
side  acceleration  (tail) 
airspeed  (coarse) 
airspeed  (fine) 
altitude 

elevator  position 

rudder  position 
yaw  rate 
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All  of  Che  calibrations  mentioned  above  were  static  in  nature,  except 
that  the  angular  accelerometer  was  exposed  to  a sinusoidal  input  because  no 
practical  method  of  providing  steady-state  angular  acceleration  was  available. 

A byproduct  of  this  calibration  method  was  rhe  capability  of  measuring  phase 
shift.  It  wa.8  assumed  that  the  oscillation  table  Itself  contributed  no  lag, 
so  that  all  phase  shift  from  input  command  to  accelerometer  signal  can  be 
attributed  to  the  accelerometer.  The  results  are  shown  in  Figure  3.  The 
transducer  apparently  develops  significant  phase  lags  at  high  frequencies, 
but  the  phase  shift  is  only  10*  to  15*  in  the  highest  frequency  region  of 
Interest.  The  phase  lead  at  very  low  frequencies  is  difficult  to  explain. 

Either  the  indication  of  phase  lead  is  an  error,  or  the  transfer  function  for 
the  transducer  is  considerably  more  complex  than  the  simple  lag  usually 
assumed.  If  the  data  found  in  these  tests  are  accurate,  it  must  be  concluded 
that  the  unit  has  been  designed  well,  since  the  phase  shift  is  zero  in  the 
frequency  region  of  greatest  interest. 

One  rate  gyro  was  selected  at  random  to  be  tested  for  phase  ohlft  in  the 
manner  described  above.  Results  are  shown  in  Figure  4.  Again,  zero  phase 
shift  seems  to  occur  in  the  frequency  range  of  most-interest,  and  a phase 
lead  of  as  much  as  10*  exists  at  low  frequencies.  Tlius  phase  shift  on  rate 
gyros  should  not  cause  any  serious  problems. 

FLIGHT  DATA  CHARACTERISTICS 

Flight  data  were  collected  for  a wide  variety  of  pilot  control  Inputs  for 
eight  flight  conditions.  For  the  purposes  of  the  current  parameter  identification 
effort,  however,  17  sets  of  time  histories,  including  only  the  two  most 
significant  flight  conditions,  were  selected.  A list  of  the  characteristics 
of  each  time  history  set  (referred  to  by  "Data  Run  Number")  is  given  in  Table  VI. 

A detailed  discussion  of  elevator  control  Inputs  is  given  elsewhere  in  this 
report.  Furthermore,  the  pilot' c record  of  the  circumstances  of  each  flight 
(from  which  data  were  taken)  are  presented  in  Table  VII.  The  details  of  the 
two  flight  conditions  are  listed  in  Table  VIII. 

The  trim  angle  of  attack  was  calculated  by  first  calibrating  the  vane 
relative  to  the  boom,  and  then  measuring  the  angle  between  the  noseboom  and 
the  wing  chord  line  (in  the  vicinity  of  the  mean  aerodynamic  chord) . Trim 
angles  of  attack  referenced  to  wing  chord  line  cannot  have  negative  values 
for  the  T-2,  yet  negative  values  are  shown  in  Table  VI.  It  seems,  therefore, 
that  an  error  exists  in  the  angle  of  attach  reference.  Such  an  error  should 
not  affect  the  perturbation  angle  of  attack  measurement  significantly. 

PILOT  CONTROL  INPUTS 

It  is  generally  recognized  (see,  for  example,  reference  (j))  that  system 
identification  is  impossible  if  the  test  input  does  not  excite  the  principal 
natural  modes  of  system  response.  Moreover,  several  investigators  (see 
references  (k)  and  (1))  have  determined  that  an  input  can  be  found  which  will 
improve  identification  performance  substantially  relative  to  traditional 
flight  test  Inputs. 
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TABLE  VII 

PILOT'S  RECORD  OF  FLIGHT  CIRCUMSTANCES 
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TABLE  VIII 

SELECTED  YT-2B  FLIGHT  CONDITIONS 


Flight  Condition  Number 

1 

2 

Nominal  Mach  Number 

0.212 

0.63 

Nominal  True  Airspeed  (ft/sec) 

236 

679 

Nominal  Pressure  Altitude  (ft) 

0 

10,000 

Gear  Position 

down 

up 

Flap  Position  (deg) 

16 

0 

Speed  Brake  Position 

closed 

closed 

Approximate  Weight  (lb) 

11,000 

11,000 

Approximate  CG  Position  (%  c) 

20 

20 

Approximate  1^  (slug-ft^) 

14,600 ' 

14,600 

Approximate  Ix  (slug-ft^) 

9,000 

9,000 

Approximate  I (slug-ft^) 
z 

19,000 

19,000 

Estimated  Trim  Angle-of -Attack  (deg) 

4.7 

1.2 

NOTE:  Weight  and  moment  of  inertia  approximations  were  based  on  reference  (i) . 


To  maximize  the  identif lability  of  the  aircraft  from  the  current  flight 
data  and  to  provide  a basis  for  the  determination  of  the  effects  of  input  type 
on  identification  performance,  it  was  decided  that  a variety  of  pilot  inputs 
would  be  used.  The  classic  step,  pulse,  and  doublet  inputs  were  included. 

The  frequencies  of  the  principal  longitudinal  modes  of  motion  (short  period 
and  phugoid)  were  estimated,  and  sine  W3ves  at  those  frequencies  were  used 
also.  The  number  of  cycles  of  the  sine  waves  was  also  varied. 

Systems  Control,  Inc.,  calculated  inputs  which  they  believed  to  be 
optimum  for  identification  of  parameters  which  have  significant  effect  on 
the  short  period  mode.  These  inputs  are  shown  in  Figures  5 and  6.  The 
shape  of  such  Inputs  is  dependent  upon  the  flight  condition  and  the 
estimated  short  period  characteristics.  On  the  theory  that  all  possible 
frequencies  should  be  represented  in  the  input,  a pseudo-white  noise  input 
was  also  selected. 
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The  actual  shape  of  the  pilot  vjput  was,  of  course,  not  always  identical 
to  the  requested  comoand  since  no  automatic  equipment  was  available  to  aid 
the  pilot  in  making  these  relatively  complex  stick  movements.  After  some 
practice,  however,  the  pilots  (who  were  not  test  pilots)  developed  a.  high 
degree  of  skill  in  executing  the  planned  inputs.  Whenever  the  real-time 
telemetry  system  was  in  operation,  the  pilot  inputs  could  be  observed  at  the 
ground  station  in  graphical  form  and  corrections  could  be  dictated  to  the 
pilot  immediately.  This  iterative  process  would  quickly  converge  to  an 
acceptable  approximation  of  the  requested  input. 

The  amplitudes  of  the  Inputs  were  also  determined  Iteratively.  Acceptable 
amplitudes  were  difficult  to  obtain  in  some  cases.  The  goal  in  the  selection 
of  input  size  was  to  force  all  significant  output  variables  to  be  large 
enough  for  a reasonable  signal-to-nolae  ratio  while  remaining  small  enough 
to  permit  the  use  of  the  linear  approximations  to  the  equations  of  motion 
aiid  to  remain  within  the  ranges  of  the  transducers.  An  additional  problem 
existed  at  flight  condition  #1,  since  the  pilots  were  reluctant  to  permit 
significant  altitude  variations  vhile  in  proximity  to  the  ground.  Unfortunately, 
altitude  variations  are  Inherent  in  the  phugold  mode,  and  it  was  therefore 
difficult  to  find  any  input  which  would  generate  a phugoid  output  without 
exceeding  any  of  the  constraints  placed  on  the  tests. 

The  rapid  doublet  input  was  eventually  abandoned  since  too  large  an 
elevator  deflection  was  required  for  sufficient  excitation  of  aircraft 
response.  The  SCI  optimum  input  was  never  successfully  generated  in  flight 
condition  #1,  ind  cannot  be  evaluated. 

SOURCES  OF  NOISE 

The  scope  of  this  effort  did  not  permit  any  attempt  at  quantitative 
analysis  of  noise  contributions  from  various  sources,  nor  of  the  sensitivity 
of  identification  performance  to  levels  of  noise  of  different  types.  It  is 
possible,  however,  to  list  the  known  sources  of  noise  and  evaluate  their 
importance  subjectively  based  on  the  experiences  of  the  investigators. 

It  is  generally  recognized  that  noise  can  be  broadly  divided  into  two 
types;  process  noise  and  measurement  noise.  In  theory,  the  modified 
Newton-Raphson  identification  technique  used  in  this  in-house  Investigation 
is  more  affected  by  process  uoise.  Given  sufficient  quantities  of  data,  this 
algorithm  will  produce  unbiased  estimates  even  in  the  presence  of  measurement 
noise,  but  the  technique  does  not  tolerate  significant  process  noise.  The 
theoretical  sensitivity  of  various  types  of  identification  techniques  to 
different  noise  types  has  been  discussed  many  times  (see,  for  Instance, 
reference  (m)),  and  will  not  be  discussed  in  this  report. 

Process  noise  can  be  further  subdivided  into  turbulence  and  modeling 
error.  Atmospheric  turbulence  was  not  noticeable  on  any  flight  at  flight 
condition  #2,  but  was  present  for  the  low  altitude  flight  condition  #1. 

Every  effort  was  made  to  fly  early  in  the  day  to  minimize  turbulence.  As 
can  be  seen  in  Table  VII,  the  flight  of  3/21/73  occurred  rather  early  in  the 
day  and  the  turbulence  level  was  such  that  the  pilot  termed  it  "none",  given 
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a choice  of  "none”,  "light",  "moderate",  and  "heavy".  On  3/14/73  the  flight 
did  not  occur  until  late  morning  and  turbulence  was  "light".  Several  additional 
flights  were  launched  at  0700  or  earlier,  but  data  from  these  flights  was 
found  unacceptable  for  reasons  other  than  turbulence  considerations.  Turbulence, 
then,  should  not  be  a factor  in  analysis  of  flight  condition  #2  data,  but 
might  have  a degrading  influence  on  identification  at  condition  #1. 

The  mathematical  model  (equations  of  motion)  used  in  this  analysis  is 
discussed  elsewhere  in  this  report.  Substantial  errors  in  the  longitudinal 
equations  per  se  are  unlikely,  but  one  must  consider  the  validity  of  the 
assumption  that  the  longitudinal  and  lateral-directional  motions  are  independent. 
Although  the  lateral-directional  motions  occurring  during  the  longitudinal 
maneuvering  are  not  large,  there  are  noticeable  excursions  in  sideslip.  It 
would  have  been  possible  to  include  in  the  model  some  sideslip  effects  and 
consider  the  sideslip  as  an  external  input  (i.e.,  an  additional  control),  but 
this  work  was  judged  to  be  outside  the  scope  of  the  current  program.  In 
summary,  the  investigators  suspect  that  process  noise  is  present,  but  not  in 
sufficient  degree  to  impede  the  identification  process,  in  spite  of  the 
theoretical  sensitivity  of  the  identification  technique  to  process  noise. 

In  theory,  every  change  in  the  nature  or  the  format  of  the  flight  motion 
data  introduces  some  measurement  noise.  It  is  clearly  desirable  to  minimize 
the  number  of  transformations  in  the  process  of  changing  phyaical  aircraft 
motions  into  digital  computer  identification  program  inputs.  However,  it 
was  necessary  that  already  available  NAVAIRDEVCEN  equipment  be  used  in  this 
investigation,  and  the  number  of  data  transformations  became  rather  large. 

The  recording  and  data  reduction  procedures  have  been  described  in  a previous 
section  of  this  report,  and  the  various  transformations  will  only  be  summarized 
here.  The  physical  motions  were  sensed  by  transducers,  which  output  electrical 
signals  in  analog  form.  These  signals  were  frequency  modulsted  and  multiplexed 
for  single  channel  recording.  They  were  then  demultiplexed,  demodulated, 
sampled  (at  20  samples/second),  and.  recorded  as  digital  information.  These 
digital  data  were  read  from  tape,  re-organized  to  a form  compatible  with  the 
NAVA IRDEV CEN  CDC  6600  computer,  decallbrated  (converted  to  physical  units), 
transformed  to  the  stability  axis  system,  and  written  onto  a disk  file  to  be 
available  for  input  to  identification  programs.  Based  on  the  investigators' 
observations  of  the  data  throughout  the  processes  outlined  above,  it  seems 
that  substantial  noise  was  Introduced  by  these  repeated  transformations. 

Specifically,  "bad  points"  arose  in  the  data  both  from  the  commutator  and 
the  analog/digital  conversion.  These  saturated  data  points  had  to  be  deleted, 
and  their  places  occupied  by  interpolated  information.  Zero  and  full-scale 
reference  values  for  both  modulation  and  analog/digital  conversion  were 
subject  to  drifting.  Each  recording/playback  process  introduced  noise  from 
the  recording  device. 

Although  no  transducer  is  perfect,  only  three  of  the  longitudinal  flight 
motion  transducers  in  the  YT-2B  caused  noticeable  problems:  elevator  position, 

airspeed,  and  angle  of  attack.  Every  transducer  selection  involves  some 
coarse/fine  trade-off.  In  the  case  of  elevator  position,  it  was  determined 
that  the  entire  range  of  deflection  would  have  to  be  included,  so  a potentiometer 
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was  Installed  on  the  actuator  and  calibrated  for  the  entire  -27°  to  +15* 
deflection  range.  At  high  dynamic  pressure  flight  conditions,  however,  the 
elevator  movement  about  its  trim  position  was  quite  small  (on  the  order  of 
magnitude  of  1°),  and  the  resulting  aignal-to-noise  ratio  for  the  data  was 
rather  poor.  No  straightforward  solution  was  found  to  this  problem  because 
the  elevator  trim  position  varied  considerably  with  flight  condition.  The 
noise  level  in  elevator  data  was  high  even  for  flight  conditions  in  which 
large  elevAtor  deflections  were  used,  indicating  that  sources  of  noise  other 
ther  resolution  difficulties  are  also  contributing. 

Airspeed  measurement  presented  a problem  analagous  to  the  elevator  position 
measurement  problem:  small  perturbations  about  a large  trim  value.  For 

airspeed  measurement,  however,  two  sensors  were  used  - a 0-500  kt  coarse 
transducer  and  a 0-50  kt  resetting  fine  transducer.  This  system  worked  well 
except  that  the  fine  transducer  was  "undecided"  about  resetting  between  48 
and  50  knots.  At,  and  just  below,  any  multiple  of  50  knots  airspeed,  the 
airspeed  measurement  becomes  useless  as  the  transducer  continuously  switches 
between  full  scale  and  zero  readings.  This  problem  is  strictly  a hardware 
phenomenon,  but  no  other  hardware  was  available. 

A filter  on  the  velocity  signal  was  used  to  process  some  of  the  data. 

This  digital  filter  simply  eliminated  velocity  values  which  were  excessively 
different  from  a previous  data  point,  and  reset  the  excessive  value  to  equal 
the  previous  data  point.  Thus  the  bursts  of  sharp  peaks  in  the  signal  were 
reduced  to  constant  value  "flap  spots".  For  the  in-house  identification 
effort,  however,  this  filter  was  applied  only  to  the  worst  data  cases  because 
extensive  use  would  have  compromised  the  principle  of  minimizing  the 
engineering  interaction  with  the  details  of  the  identification  process. 

Angle  of  attack  was  measured  by  a vane  mounted  on  a noseboom.  The  vane 
mounting  position  was  sufficiently  forward  of  the  wing  for  upwash  effects  to 
be  negligible.  However,  the  vane  appeared  to  have  dynamic  characteristics 
which  varied  from  day  to  day  as  the  friction  on  the  pivot  changed  with  weather, 
time,  etc.  Unfortunately,  no  means  of  dynamic  calibration  for  the  vane  was 
available. 

In  summary,  the  investigators  suspect  that  most  of  the  noise  in  the  flight 
data  was  measurement  noise.  It  is  especially  significant  that  measurement 
noise  was  considerable  in  the  control  input,  since  theoretical  parameter 
identification  studies  generally  consider  the  input  as  a noise-free  measurement. 


SELECTION  OF  MATHEMATICAL  MODEL 

The  general  equations  of  motxon  for  aircraft  flight  are  quite  complex,  but 
numerous  assumptions  are  often  made  to  justify  simplification  of  the  equations 
for  a particular  application.  Two  simplifying  assumptions  were  implicit  in 
the  definition  of  this  project:  (1)  longitudinal  equations  may  be  decoupled 

from  lateral-directional  equations , and  (2)  small  perturbation  theory  and 
small  angle  approximations  are  applicable  to  the  flight  motions.  Equations 
based  on  these  assumptions  are  deve..oped  in  many  texts,  but  reference  (n) 
was  used  for  this  project. 
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The  modified  Newton-Raphson  computer  program  (and  moat  other  identification 
schemea)  requlrea  that  the  model  be  in  the  atate  vector  form 

x ■ Fx  + Gu 

where  x la  the  vector  of  the  motion  atatea,  u la  the  control  vector,  F la  the 
matrix  of  atabillty  derivatives,  and  G ia  the  matrix  of  control  derivatlvea. 

The  equatlona  of  motion  are  not  usually  written  in  such  a form  when  used  for 
dynaarf.cs  analysis,  but  they  can  be  manipulated  to  the  proper  form.  The 
resulting  equations  are  listed  In  Table  IX. 

During  the  identification  process,  the  equations  were  simplified  even 
further,  as  z«.  v X«e,  and  were  all  taken  to  be  zero.  The  quantities 

u and  ff  were  not  measured,  but  were  calculated  from  a and  a using  auxllliary 
equations. 

With  the  equations  of  amtion  in  the  atate  vector  form,  it  is  clear  that 
no  straightforward  aolutlon  for  is  possible.  Various  schemes  were  considered, 

including  the  use  of  the  auxllliary  relationship  M*  ■ tail  length  X Z*  and 

obtaining  data  at  large  V , to  permit  separation  of  from  the  lumped  parameters, 

but  all  were  eventually  discarded  as  grossly  inaccurate  or  impractical.  To 
separate  the  identified  derivatives  for  comparison  to  estimates,  either  the 
estimate  can  be  assuamd  correct  or  the  estlamttes  can  be  combined  into  lumped 

parameters  and  the  comparisons  done  with  those  quantities. 


DEVELOPMENT  OF  A PRIORI  PARAMETER  SET 

Despite  occassional  claims  that  various  parameter  identification  techniques 
will  function  without  any  a priori  or  start-up  parameter  estimates,  it  ia 
generally  agreed  that  the  use  of  such  estimates  substantially  increases  the 
chances  of  success.  A priori  paraamter  estimates  which  merit  a high  level  of 
confidence  may  also  serve  as  a basis  of  coeq>arlaon  for  results,  but  should 
never  be  considered  the  "right  answers".  Indeed,  if  a priori  estimates  were 
to  be  considered  the  "right  answers",  then  there  would  be  no  need  for  paraswter 
identification  from  flight  data. 

A priori  parameter  estimates  for  the  T-2  for  the  flight  conditions  of 
current  interest  were  developed  from  information  in  references  (f ) , (i) , and 
(o).  The  methods  used  by  the  aircraft  manufacturer  to  generate  this  original 
information  are  not  known,  but  a reasonable  level  of  confidence  would  seem 
appropriate  for  any  information  released  to  the  public  by  the  aircraft  design 
team. 

The  estimated  parameters  are  given  in  Table  X in  several  forms.  Only  one 
set  of  information  is  presented,  but  several  formats  were  used  as  a convenience 
to  the  user.  This  information  has  been  supplied  to  all  parties  receiving  the 
flight  data. 
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TABLE  X 

A PRIORI  VALUES  FOR  YT-ZB  AERODYNAMIC  PARAMETERS 


Parameter 

Dimensions 

Flight 

1 

Condition 

2 

c 

0.642 

0.106 

L 

C 

• 

0.110 

0.018 

D 

CLry 

- 

4.62 

5.50 

CD,-y 

- 

0.172 

0.228 

CL,f 

1/rad 

1.049 

- 

CL, 

1/rad 

0.486 

0.401 

e 

°Lu 

- 

0 

0.031 

CDu 

- 

0 

0 

C°u 

- 

0 

0.011 

m 

-0.597 

-0.429 

nVv 

S; 

- 

-3.95 

-5.50 

% 

- 

-10.59 

-11.86 

1/rad 

-1.13 

-0.994 

X 

1/sec 

-0.046 

-0.016 

u 

X 

fy 

ft/(secz  - rad) 

-23.23 

-36.91 

z 

u 

Z 

fy 

1/sec 

-0.269 

-0.122 

ft/(sec^  - rad) 

-233.75 

-1669.24 

Z*e 

ft/(sec2  - rad) 

-24.02 

-121.30 

Mu 

l/(ft  - sec) 

0 

0.0017 

**« 

l/(sec2  - rad) 

-5.11 

-22.50 

M* 

l/(sec  - rad) 

-0.53 

-1.57 

\ 

l/(sec  - rad) 

-1.42 

-3.39 

Mfie 

l/(sec2  _ rad) 

-9.68 

-52.12 

-g  (cos  V ) /U 

1/sec 

-0.134 

-0.047 

Z„  V<"o  - 

1/sec 

-0.27 

-0.12 

V<“o  - V 

1/sec 

-0.974 

-2.46 

Wo  + V/(Uo  • 

- 

1.0 

1.0 

-g  (sin  vo)/(U0 

-z;> 

1/sec 

0 

0 
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TABLE  X (CON'T) 

A PRIORI  VALUES  FOR  YT-2B  AERODYNAMIC  PARAMETERS 


Parameter 


Dimensions  Condition 


X AJ  I/sec 

(V  o 

X*e/U0  1/sec 

(uo  c-‘u + M;  ■ z;»  1/sec2 

+ Z„/W0  - Z.»  l/.ec2 

Mq  + M;  «“o  + Z,)/ff0  - Zp)  !/•« 

-(M^  g Sin  vo)/(UQ  - Z^)  1/sec2 

Z5e/(Uo  " Z*>  1/sec 

M5e  + (Zfie  M;/(Uq  - Z^)  l/sec2 


■0.097 

-0.054 

0 

0 

1.43 

1.34 

•4.59 

-18.6 

1.95 

-4.96 

0 

0 

0.102 

-0.18 

9.63 

-51.8 

43 


NA  DC-74181-30 


The  conditions  of  each  actual  data  run  are,  of  courae,  slightly  different 
from  each  other  and  from  the  nominal  flight  condition.  For  the  In -house 
Identification  analyals,  a separate  set  of  combined  start-up  parameters  was 
generated  for  each  run  reflecting  the  differences  In  trim  velocity. 


IN-HOUSE  IDENTIFICATION 
COMPUTER  PROGRAM 


ACQUISITION 

The  Identification  technique  used  in  this  effort  was  a modified  Newton - 
Raphson  algorithm.  The  computer  program  was  provided  (Informally)  by  Mr.  Larry 
Taylor  of  NASA-Langley  Research  Center.  The  rationale  of  the  program  and  an 
example  of  its  application  are  given  In  reference  (p).  Reference  (q)  Is  a 
"user's  guide"  for  the  program.  In  fact,  the  computer  program  has  been 
documented  so  thoroughly  that  It  would  be  redundant  to  discuss  It  in  this 
report.  Minor  changes  were  made  to  the  program  In  accordance  with  references 
(r)  and  (a),  and  some  modifications  had  to  be  made  to  adapt  the  program  to  the 
current  NAVAIRDEVCEN  CDC  6600  cosqmter  operating  system. 

VERIFICATION 

At  the  time  of  its  acquisition  by  NAVAIRDEVCEN,  the  computer  program  had 
been  configured  for  analysis  of  lateral-directional  dynamics.  To  minimize 
the  number  of  immediate  programing  changes,  it  was  decided  that  the  program 
would  first  be  applied  to  lateral-directional  "flight  data"  generated  by  a 
digital  computer  simulation  of  the  T-2,  and  then  the  program  would  be  tried  on 
longitudinal  data.  Results  of  the  first  run  on  lateral-directional  data 
are  shown  in  Table  XI.  Since  the  data  are  theoretically  noiseless  and  the  same 
equations  of  motion  were  used  In  both  the  simulation  and  identification 
computer  programs,  the  identification  should  be  perfect  if  all  controls  are 
used.  Both  aileron  and  rudder  were  applied  in  the  simulation  program. 

Run  1 resulta  indicate  that  the  parameter  identification  was  excellent, 
but  not  perfect.  Small  errors  existed  in  almost  all  parasieters,  and  N^,  Yp, 

and  were  not  identified  at  all.  It  is  customary  to  neglect  Y*  and  Y 

In  most  applications,  and  N^  was  very  small  In  this  case,  although  Np  can  be 

an  important  derivative  in  some  cases.  It  Is  clear,  then,  that  some  noise  was 
present  after  all.  The  investigators  suspect  that  the  source  of  the  noise 
was  an  incoaipatibility  in  the  Integration  schemes  used  in  the  digital  computer 
programs  for  identification  and  flight  simulation.  Any  quantification  error 
resulting  from  integration  algorithm  mismatch  was  a phenomenon  unique  to 
application  of  the  identification  method  to  digitally  simulated  flight  data, 
so  that  the  problem  should  not  exist  in  applications  of  the  program  to  true 
flight  data,  and  no  further  investigation  waa  made. 

Several  types  of  gusts  were  added  to  the  flight  data  to  test  the  reaction 
of  the  identification  program  to  process  noise.  The  characteristics  of  runs 
2 through  5 are  given  in  Table  XII . It  can  be  seen  from  the  results,  given 
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TABLE  XII 

NOISE  CHARACTERISTICS  OF  COMPUTER-GENERATED 
T-2  FLIGHT  TIKE  HISTORIES 


Run 

No. 


Characteristics 


1.  Noiseless.  Should  be  perfect  identification.  Used  as  reference. 

2.  Low  amplitude,  low  frequency,  sinusoidal  g and  r gusts.  No  p gusts. 

3.  Moderate  amplitude,  moderate  frequency,  sinusoidal  0 and  r gusts; 
and  low  amplitude,  high  frequency  p gusts. 

4.  Amplitude  and  spectral  characteristics  of  8 gusts  per  MIL-F-8785B. 

No  r or  p gusts. 

5.  Amplitude  and  spectral  characteristics  of  8 and  r gusts  per 
MIL-F-8785B.  No  p gusts. 

6.  No  gusts.  Lag  with  .098  second  time  constant  Introduced  in  3 vane. 

7.  No  gusts.  Pseudo-random  noise  in  all  measurements.  Noise  amplitudes 
selected  on  basis  of  Cornell  Aeronautical  Laboratory  experience  in 
X-22  identification  program. 
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in  Table  XI,  chat  the  the  quality  of  tho  identification  process  deteriorates 
rapidly  as  gust  intensity  is  increased.  The  presence  of  turbulence,  having 
the  characteristics  set  forth  in  section  3.7  of  military  specification 
MIL-F-8785B  (ASG)  for  ambient  atsiospheric  turbulence,  is  sufficient  to 
seriously  degrade  the  performance  of  the  identification  method.  The  results 
obtained  by  the  computer  program  from  flight  data  acquired  under  such  turbulence 
conditions  could  be  used  only  to  specify  ±10%  estimates  of  Lg,  Lfta,  L6r,  N0, 

and  Njr,  a ±15°  estimate  of  L^,  and  no  estimates  of  the  other  parameters. 

These  poor  results  are  not  unexpected,  since  the  modified  Newton-Raphson 
algorithm  is  known  to  give  biased  results  in  the  presence  of  process  noise. 
Whenever  a consistent  set  of  turbulence  characteristics  of  a particular 
level  (e.g.,  light  turbulence,  moderate  turbulence,  MIL-F-8785  ambient 
turbulence)  is  present,  the  p gusts  are  far  more  damaging  to  identification 
quality  than  r or  p gusts. 


Run  6 was  made  to  evaluate  the  effect  of  a lag  in  the  p-vane,  and  all 
other  noise  was  testporarily  deleted  to  Isolate  this  effect.  The  presence  of 
this  small  uncompensated  time  lag  in  the  sideslip  vane  seriously  corrupts 
estimation  of  Yg,  Yftr,  Nf,  and  Lf,  and  has  a lesser  effect  on  the  other 

derivatives.  This  result  is  curious,  since  the  four  listed  derivatives  are 
not  those  which  one  would  expect  to  be  affected  most,  but  the  important 
conclusion  is  that  the  existence  of  the  lag  is  sufficient  to  seriously 
disrupt  the  identification  process. 


In  theory,  the  operation  of  the  modified  Newton-Raphson  algorithm  is  not 
biased  by  the  presence  of  zero  mean  white  Gaussian  measurement  noise.  While 
the  results  from  the  application  of  the  identification  method  to  any  individual 
set  of  time  histories  may  be  corrupted  by  the  presence  of  measurement  noise, 
the  mean  value  of  the  parameter  estimates  will  approach  the  true  parameter 
values  during  repeated  analysis  of  different  sets  of  time  histories  for  the 
same  flight  condition.  No  source  of  truly  white  noise  was  available  for  data 
generation,  so  a zero-mean  approximation  was  developed.  The  noise  amplitudes 
were  selected  on  the  basis  of  Calspan  Corporation  experience  with  flight 
testing  of  the  X-22  aircraft  in  reference  (e) . 


The  results  from  the  identification  program  operating  in  the  presence 
of  such  pseudo-random  measurement  noise  may  be  described  as  fair.  The 
Lp,  Lg,  L6a,  L^r,  Nf,  Ng,  N$a,  and  parameters  are  estimated  within  ±10% 

error.  These  results  sre  presented  as  run  7 in  Table  XI. 

The  longitudinal  form  of  the  identification  program  was  also  tested, 
albeit  briefly.  When  the  initial  parameter  values  are  identical  to  the 
correct  results,  the  first  iteration  of  the  program  has  a very  low  fit  error. 

If  the  program  is  allowed  to  continue  to  iterate,  a slightly  smaller  fit  error 
is  obtained  with  a set  of  stability  derivatives  somewhat  different  from  the 
correct  value.  Presumably,  this  problem  is  analogous  to  the  slight  disparities 
found  in  the  lateral-directional  data  identification. 
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It  oust  be  noted  that  the  weighting  matrix  in  the  identification  program 
was  not  used  during  these  applications  of  the  program.  These  experiments  were 
performed  immediately  after  acquisition  of  the  program,  and  it  was  felt  at 
that  time  that  use  of  the  weighting  matrix  would  be  difficult  and  would 
violate  the  concept  of  an  "off-the-shelf"  program  to  be  used  in  "cookbook" 
fashion.  Later  experiences,  described  elsewhere  in  this  report,  revealed 
that  the  weighting  matrix  could  be  used  easily  and  in  straightforward  manner, 
and  that  its  use  would  substantially  improve  the  performance  of  the  program 
on  noisy  data, 

IN-HOUSE  IDENTIFICATION  PROCEDURES 

A number  of  "data  runs"  were  selected  from  the  original  flight  data. 

These  runs  were  processed  by  program  MASSAGE  (see  listing.  Appendix  A)  to 
a reduced  form.  The  runs  were  then  analyzed  individually  by  program  NEWTON 
(see  listing.  Appendix  B).  In  most  cases  3S0  data  points  were  extracted  from 
the  input  data  file  for  each  identification  run.  The  time  location  of  these 
points  was  selected  to  place  the  control  input  several  seconds  from  the 
beginning.  For  some  runs,  less  than  350  points  were  available  (data  run  6. 
for  Instance,  only  had  200  data  points  before  occurrence  of  excessive  noise). 

At  20  data  points  per  second,  these  350  point  data  sets  were  17.5  seconds 
in  length.  In  an  attempt  to  extract  velocity  derivatives  from  longer  data 
runs,  a mechanism  was  provided  in  the  program  for  reading  only  every  second, 
or  third,  or  fourth  data  point,  so  that  350  data  points  were  still  used,  but 
the  time  interval  was  multiplied. 

A bad  point  filter  was  used  in  the  identification  program,  which  substituted 
a reasonable  value  (calculated  by  extrapolation  from  previous  points)  for  any 
data  point  which  was  completely  off  scale.  A somewhat  more  restrictive  filter 
was  applied  to  airspeed  in  those  instances  where  spikes  existed  from  transducer 
switching.  This  filter  also  operated  with  the  substitution/extrapolation 
method.  These  were  the  only  filters  used  on  the  flight  condition  #1  data  - 
there  was  no  general,  lowpass  filtering  done  during  A/D  conversion. 

The  very  low  slgnal/nolse  ratios  in  the  elevator  flight  data  for  flight 
condition  #2  led  to  the  use  of  several  relatively  extreme  measures.  The 
elevator  signal  was  treated  with  a weighted  average  smoothing  routine.  Only 
a three  point  width  smoothing  interval  was  used  to  prevent  excessive  loss  of 
high  frequency  signal  content,  and  the  signal  was  smoothed  either  once  or  twice. 
The  suppression  of  the  elevator  signal  to  zero  following  the  significant 
portion  of  the  control  input  was  also  done  for  the  data  sets  containing 
"spikes".  The  smoothing  was  done  after  the  suppression  to  maintain  a smooth 
transition. 

The  identification  program  contained  a weighting  matrix  which  weighted 
the  importance  of  the  time  history  matches  in  the  parameter  update  iteration 
process.  The  program  printout  lists  the  weighted  contributions  to  the  overall 
time  history  fit  error,  and  the  diagonal  elements  of  the  weighting  matrix 
should  each  be  the  inverse  of  the  corresponding  weighted  fit  error  term.  This 
process  of  adjusting  the  weighting  matrix  thus  became  an  outer  iterative  loop  - 
the  inner  loop  being  the  parameter  determination  with  a fixed  weighting  matrix. 
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The  weighting  matrix  was  updated  manually,  and  after  two  or  three  iterations 
the  weighted  fit  error  terms  would  usually  all  be  nearly  equal  to  1 (between 
0.9  and  1.1),  which  of  course  was  a terminal  condition  since  further  inversion 
would  yield  no  changes.  The  process  of  correcting  the  weighting  matrix 
clearly  could  be  Included  in  the  program  itself  if  one  wished  to  make  the 
necessary  effort.  The  effect  of  the  weighting  matrix  was  to  compensate  for 
the  dimensional  differences  in  the  time  histories  (in  effect  a normalization) 
and  to  proportion  the  significance  of  each  time  history  inversely  with  its 
noise  level. 

In  most  cases,  all  the  Identifiable  derivatives  were  left  as  variables, 
but  on  the  longer  time  (such  as  52.5  seconds)  runs,  the  short  period 
derivatives  were  held  constant  at  the  values  determined  for  those  derivatives 
by  analysis  of  the  .cirst  17.5  seconds  of  the  same  data  run. 

The  variables  read  into  the  identification  program  from  the  data  massage 
program  were  A , r*,  p,  q,  q (calculated  from  linear  accelerations),  (measured 

directly),  u,  a , «z,  h,  *b>  Of  these  variables,  q,  h,  and  *b  were  not 

usually  used  by  the  identification  program.  The  q signal  appeared  far  more 
noisy  than  the  q^  signal  and  therefore  was  abandoned.  The  mechanization  of 
the  computer  program  had  no  provision  for  use  of  altitude  information,  and 
speed  brakes  simply  were  not  used  by  the  pilot  during  any  of  the  selected 
data  runs. 

Since  the  equations  of  motion  in  the  program  are  in  state  variable  format, 

a and  a were  used  to  calculate,  respectively,  u and  *>,  which  were  the  matched 
X z 

variables  along  with  o,  u,  8,  q,  and  q. 

According  to  reference  (q) , the  output  from  the  computer  program  includes 
"the  approximate  standard  deviation  if  the  inverse  of  the  noise  covariance 
matrix  is  used"  for  the  weighting  matrix.  The  investigators  on  this  project 
felt,  however,  that  the  method  of  calculation  of  this  "approximate  standard 
deviation"  was  not  necessarily  valid.  It  was  decided  not  to  attempt  to  use 
this  quantity  as  a measure  of  estimation  quality,  and  no  other  confidence 
level  Indicator  was  developed. 

IN-HOUSE  IDENTIFICATION  RESULTS 
SUMMARY  OF  FLIGHT  CONDITION  NO.  1 

Flight  condition  #1  time  history  data  included  data  runs  2 through  6.  A 
tabular  summary  of  the  stability  derivatives  obtained  from  the  various 
modified  Newton-Raphson  identification  attempts  with  these  data  is  given  in 
Appendix  C.  A tabular  summary  of  the  results  of  significant  runs,  reduced  with 
respect  to  velocity  and  reduced  with  an  assumed  and  Z^,  is  given  in 

Table  XIII.  The  reduced  form  will  be  used  for  most  discussions  because  it  is 
more  universal  and  more  easily  compared  to  estimates,  even  though  some 
inaccuracy  is  clearly  introduced  by  the  assumption  of  a known  M..  Although 
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the  trim  velocity  has  been  eliminated  as  a direct  effect  on  the  identified 
derivatives,  it  is  still  true  that  velocity  varies  significantly  among  the 
five  data  runs,  and  some  variation  in  the  dimensional  stability  derivatives 
is  to  be  expected  with  variation  in  trim  velocity.  Time  history  matches  for 
all  of  the  runs  included  in  Table  XIII,  as  well  as  some  other  matches  from 
preliminary  runs,  are  shown  in  Figures  7 through  16. 

EFFECT  OF  WEIGHTING  MATRIX  AND  START-UP  VALUES 

The  first  identification  attempt  was  made  on  data  run  3 with  an  identity 
weighting  matrix  and  all  zero  start-up  values,  but  only  a divergence  was 
obtained.  A converging  run  was  eventually  obtained  using  a set  of  start-up 
values  which  were  the  "best  guess"  a priori  stability  derivatives,  based  on 
the  trim  velocity  of  data  run  2 (which  had  also  been  tried  with  zero  start-up 
values  but  diverged).  Under  these  circumstances,  the  program  converged,  but 
the  stability  derivatives  were  virtually  unchanged  from  the  start-up  values 
and  the  time  history  matches  (see  Figure  7)  showed  significant  problems.  Also 
visible  in  Figure  7 are  the  poor  slgnal-to-noise  ratio  for  the  elevator 
input  and  a typical  unfiltered  velocity  transducer  switching  problem.  This 
first  identification  run,  however,  did  provide  sufficient  information  to 
calculate  a weighting  matrix.  The  next  run  (Figure  8)  showed  significant 
improvement,  and  a second  iteration  (Figure  9)  brought  about  a respectable 
time  history  match  and  a series  of  weighted  fit  error  contributions  nearly 
equal  to  one. 

EFFECT  OF  ZEROED  CONTROL  INPUT 

In  an  effort  to  assess  the  damage  caused  to  the  identification  process  by 
the  noisy  elevator  input,  a run  was  made  with  Ag  set  to  zero  after  the  pii'-  . 

control  input  was  clearly  completed.  Unfortunately,  an  overzealous  velocity 
filter  set  questionable  velocity  data  to  zero  as  well,  rather  than  simply 
smoothing.  This  distortion  of  the  velocity  time  history  was  felt  to  be 
relatively  insignificant,  however.  The  time  history  matches  (see  Figure  10) 
are  quite  similar,  but  comparison  of  run  DR3-1W2  and  run  DR3-1WDEZ  results  in 
Table  XIII  shows  that  noticeable  changes  occur  in  the  identified  derivatives, 
so  control  noise  clearly  was  a factor. 

One  identification  run  was  accidentally  made  with  the  control  input  zeroed 
for  the  entire  time  history  with  the  time  histories  of  ihe  other  variables 
intact.  Such  an  exercise  is  only  of  academic  interest,  of  course,  but  the 
results  do  reveal  some  interesting  aspects  of  the  operation  of  the  identification 
algorithm.  The  program  returned  values  for  the  control  derivatives  which  were 
equal  to  the  a priori  values,  rather  than  calculating  Very  large  numbers  or 
refusing  to  find  a solution.  The  implication  is  clearly  that  the  algorithm  tends 
to  retain  the  a priori  parameter  values  when  given  insufficient  information 
to  estimate  the  parameters.  This  behavior  may  explain  the  phenomenon!  of 
reasonable  estimation  of  the  Gpeed  derivatives  from  short  data  records  with 
little  or  no  speed  variation. 
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KEY  TO  TIME  HISTORY  COMPARISONS 

Solid  lines  are  flight  data.  Dotted  lines  are  calculated  responses.  For 
elevator  time  history,  solid  line  is  measured  input  data  while  dotted  line 
is  smoothed  version  actually  used  as  input  for  identification  program. 

DE  - elevator  deflection,  degrees,  positive  trailing  edge  up 

A - angle  of  attack,  degrees,  positive  for  airplane  nose  up 

U * airspeed,  ft/sec 

TH  - pitch  angle,  degrees,  positive  for  airplane  nose  up 

Q ■ pitch  rate,  rad Ians /seconds,  positive  clockwise  viewed  from  port  side 
of  aircraft 

AX  ■ longitudinal  acceleration,  ft/sec^,  positive  forward 
AZ  ■ normal  acceleration,  ft/sec^,  positive  down 

QD  - pitch  acceleration,  radlans/sec^,  positive  clockwise  viewed  from  port 
side  of  aircraft 

TIME  ■ elapsed  time,  seconds 

Run  number  coding: 

1.  Digit  after  "DR"  is  data  run  number. 

2.  "U"  after  hyphen  indicates  that  only  speed  derivatives  were  variables. 

3.  First  digit  is  time  length  in  multiples  of  17.5  seconds. 

4.  Digit  after  "W"  is  number  of  iterations  to  obtain  satisfactory  weighting 
matrix. 

5.  Digit  after  "S"  is  number  of  times  smoothing  applied  to  input  data. 

6.  Absence  of  "W"  or  "S"  indicates  absence  of  weighting  or  smoothing, 
respectively. 

7.  "DEZ"  indicates  elevator  set  to  zero  following  significant  control  input. 

■j 
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FIGURE  7.  TIME  HISTORY  COMPARISON  DR3-1 


FIGURE  9.  TIME  HISTORY  CGMHUUSON  DR3-1W2 


FIGURE  10.  TIME  HISTORY  COMPARISON  DR3-1W1DEZ 


FIGURE  11.  TIME  HISTORY  CCMIARISON  DR3-2W2 
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FIGURE  12.  TIME  HISTORY  COMPARISON  DR3-U3W2 


FIGURE  14.  TIME  HISTORY  COMPARISON  DR4-1W4 
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FIGURE  14  (COOT).  TIME  HISTORY  CQMBkRXSON  DR4-1W4 
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If  the  algorithm  does  tend  to  keep  parameters  at  their  a priori  values  when 
faced  with  insufficient  data,  then  parameter  results  would  be  a strong 
function  of  their  start-up  values  if  the  flight  data  were  Incomplete  or 
inappropriate.  This  theory  was  not  tested  because  all  radical  changes  in 
start-up  values  lead  to  divergence. 

EFFECT  OF  HOLDING  SHORT-PERIOD  DERIVATIVES  CONSTANT 

Utilising  a 35  second  section  of  data  run  3,  two  different  identification 
procedures  were  used.  One  run  (DR3-2W2)  had  all  derivatives  as  variables 
while  the  other  (DR3-U2W1)  had  short  period  derivatives  held  fixed  at  the 
values  determined  by  the  best  17.5  second  run  (DR3-1W2)  and  only  the  speed 
derivatives  variable. 

Only  11^  was  significantly  different  in  the  two  sets  of  results.  As  can 

be  seen  in  Table  XIII,  the  speed  derivatives  were  essentially  identical.  The 
time  history  match  (Figure  11)  is  for  run  DR3-2W2,  but  an  identical  plot  was 
generated  for  DR3-U2W1.  Both  time  history  matches  were  disturbed  by  apparent 
atmospheric  turbulence.  At  approximately  10  seconds  and  28  seconds  significant 
variations  in  angle  of  attack  occur  without  any  visible  change  in  * , Pitch 

rate  and  pitch  angle  also  change  at  this  time,  but  not  as  severely  as  angle  of 
attack.  On  the  52.5  second  run  with  the  same  data  (see  Figure  12)  even  more 
turbulence  is  obvious,  and  a large  disturbance  in  velocity  is  visible  lagging 
significantly  behind  angle  of  attack.  The  model  equations  of  motion  have 
neither  gust  inputs  nor  nonlinear  terms,  so  the  model  insists  on  creating 
smooth  time  histories.  The  velocity  time  history,  for  Instance,  rejoins  the 
flight  data  after  the  gust  disturbance. 

DATA  RUN  #3 

Several  identification  attempts  using  data  run  3 have  already  been 
discussed  for  other  purposes,  but  it  remains  necessary  to  take  an  overall 
look  at  data  run  3 results.  The  principal  results  (DR3-1W2)  agree  well  with 
the  a priori  values  except  for  X^,  which  is  negative  but  much  smaller,  and 

M*e,  which  is  about  251  below  the  a priori  value.  The  other  runs  (DR3-1WDEZ 

and  DR3-2W2)  which  had  variable  short  period  derivatives  had  a similarly  low 
Mfte,  and  higher  than  the  a priori  value  (and  about  25 1 higher  than  DR3-1W2) , 

and  a widely  varying  X . The  speed  derivatives  are  surprisingly  consistent 

or 

from  the  short  runs  to  the  long  runs,  and  also  similar  to  the  a priori  values. 

The  time  history  match  (Figure  9)  for  angle  of  attack  seems  to  be  both 
slightly  out  of  phase  and  offset  somewhat.  The  offset  may  be  due  to  a poorly 
defined  trim  condition  on  angle  of  attack  due  to  the  atmospheric  turbulence. 

The  time  history  matches  for  the  other  variables  are  quite  good  except  for  the 
previously  discussed  areas  of  turbulence. 
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DATA  RUN  n 

As  can  be  seen  from  Table  XIII , the  results  from  data  run  2 are  significantly 
different  from  the  other  data  runs.  Specifically,  Zy  has  an  incorrect  sign 

(although  it  would  be  foolish  to  expect  to  estimate  speed  derivatives  from  a 
short  data  record),  and  are  unusually  large,  and  M/\e,  although  smaller 

than  the  a priori  value,  is  much  larger  than  that  obtained  from  any  other  runs. 
Run  2 has  a higher  trim  velocity  than  any  other  data  run,  but  this  factor  alone 
is  not  enough  to  explain  the  difference.  Considering  C^  as  Independent  of 

trim  velocity,  for  instance,  Z for  data  run  2 should  be  -275.6  (based  on  a 

value  of  -235.3  for  Z from  data  run  3)  but  was  estimated  at  -311.8. 

Figure  13  shows  the  time  history  matches  for  data  run  2,  including  ax, 

az,  and  q in  addition  to  the  »>t  u,  4,  and  q shown  for  data  run  3.  The 

angle  of  attack  time  history  match  is  far  from  perfect.  Some  phase  shift 
seems  to  exist,  with  the  modeling  leading  the  flight  data,  and  the  model 
response  is  smaller  in  amplitude  than  the  flight  data.  The  flight  data  is 
nearly  symmetric  with  respect  to  positive  and  negative  excursions  in  angle  of 
attack;  the  model  response  has  smaller  negative  excursions  than  positive, 
creating  a noticeable  fit  error. 

The  time  history  match  is  not  a match  at  all  - the  model  response  has 

a large  lead  relative  to  the  flight  data  and  the  principal  response  is  in 
the  opposite  direction.  As  stated  previously,  the  program  performs  a match 
on  u,  not  on  a^  itself,  and  u contains  a g°  term  which  in  this  case  is  much 

larger  in  magnitude  than  a^.  Thus  the  u match  is  largely  a repeat  of  the  8 

match  for  this  data  set.  The  opposite  sign  trend  in  the  response  nevertheless 
motivated  a thorough  investigation  of  sign  convention  relative  to  a , but 
no  mistakes  were  found. 

The  a^  situation  is  similar  to  the  a^  problem  mentioned  above,  in  that  the 

match  is  partially  a repeat  of  the  q match  rather  than  a alone,  but  the  a 

z z 

time  history  match  is  really  quite  good.  Of  course  the  level  of  excitation  is 
quite  different  - a is  approximately  ±30  ft/sec^  while  a is  ±2  ft/sec^. 

Z X 

The  q time  history  comparison  is  quite  unusual,  since  the  model  response  is 
much  noiscr  than  the  flight  data.  The  flight  data  for  q are  taken  from  the 
angular  accelerometer,  which  was  less  noisy  than  expected.  The  other  q 
measurement,  calculated  from  nose  and  tail  linear  accelerometers,  was  far  ton 
noisy  to  be  of  any  use.  The  model  response,  which  is  directly  proportional  to 
elevator  deflection,  is  quite  noisy  since  the  elevator  deflection  is  so  noisy. 
Although  the  measured  fit  error  is  quite  large  for  this  time  history  pair,  the 
trends  are  clearly  correct  and  it  seems  fair  to  judge  that  q is  matched. 
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DATA  RUN  #4 

The  stability  derivatives  extracted  from  data  run  4 are  quite  similar  to 
the  a priori  values  (and  hence  to  data  run  3 results) . The  widely  varying  X^ 

has  essentially  the  same  value  here  as  for  data  run  2.  is  unusually  small 

for  data  run  4,  however,  being  87.  below  the  a priori  value.  The  value  of 

is  noticeably  (107.)  below  the  a priori  value,  but  it  is  very  close  to  the 
results  from  previous  data  runs.  The  M*e  extracted  for  this  run  is  a compromise, 

almost  evenly  spaced  between  the  results  of  run  3 and  run  2.  The  extracted 
derivatives  are  listed  in  Table  XIII. 

Time  history  comparisons  for  data  run  4 are  shown  in  Figure  14.  The 
comments  made  for  data  run  2 (above)  apply  almost  exactly  to  data  run  4 as  well. 
Angle  of  attack  time  shift  and  amplitude,  a time  shift  and  direction  of  excur- 
sion, and  q model  response  noise  were  stillproblems. 

DATA  RUN  #5 

The  control  input  for  data  run  5 was  similar  to  that  of  data  run  2,  except 
that  the  former  was  somewhat  lower  in  frequency  and  amplitude  than  the  latter. 

The  trim  airspeed  for  data  run  5 is  11  ft/sec  less  than  data  run  2,  which  could 
lead  to  some  differences  in  parameter  values,  but  the  data  run  5 trim  airspeed 
is  essentially  in  the  center  of  the  range  of  trim  airspeeds  covered  by  the 
data,  so  overall  comparisons  of  results  without  consideration  of  airspeed 
should  be  possible. 

The  most  obvious  feature  of  the  results  for  data  run  5 is  that  both  X and 

u 

X^  are  estimated  to  be  positive,  while  all  other  estimates  for  these  derivatives 

were  negative.  The  estimate  for  M is  somewhat  large,  but  not  as  large  as 

the  estimate  from  data  run  2.  The  M*  estimate  falls  near  the  center  of  the 

e 

range  defined  by  previous  estimates. 

The  time  history  matches  for  data  run  5,  shown  in  Figure  15,  appear  to  be 
slightly  different  from  those  for  previous  runs.  The  angle  of  attack  phase  and 
amplitude  problems  are  slightly  smaller  than  for  the  runs,  but  the  a z fit  has 

deteriorated  noticeably.  The  ax  flight  data  seems  to  have  a steady-state 

offset,  but  the  model  and  flight  data  excursions  seem  to  be  approximately 
equal  in  magnitude,  identical  in  direction,  and  nearly  in  phase. 

DATA  RUN  #6 

Excessive  noise  rendered  the  latter  portion  of  data  run  6 useless,  so  only 
200  data  points  were  used.  The  results  for  the  shortened  run  6 are  given  in 
Table  XIII.  A reasonable  value  for  X was  obtained,  while  X was  exactly  zero. 

The  remainder  of  the  stability  derivatives  were  very  similar  to  the  results 
from  previous  runs.  Although  it  may  seem  unlikely  that  the  speed  derivatives 
could  be  estimated  from  a 10  second  data  run,  the  magnitude  of  the  velocity 
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perturbation  la  actually  larger  (about  20  ft/aec)  In  data  run  6 than  In 
aeveral  other  data  run a. 

Time  hlatory  comparisons  for  data  run  6 ara  shown  In  Figure  16.  The 
angle  of  attack  tlae  hlatory  match  la  In  phase,  but  a aerloua  amplitude  error 
exists.  The  model  reaponae  la  approximately  75X  of  the  ln-flight  reaponae. 
The  a match  la  good.  The  a model  reaponae  leada  the  flight  data  alight ly 

and  the  excuralon  la  once  again  In  tha  oppoalta  direction. 

INTERPRETATION  OF  X RESULTS 

a 

Two  atablllty  derivatives,  and  M^,  vary  more  widely  than  anticipated 

In  the  reaulta  dlacuaaed  above,  and  It  would  be  dealrable  to  Identify  reaaona 
for  thla  behavior.  Although  thla  report  generally  daala  In  dlmenalonal 
derlvatlvea,  which  are  conaldered  more  Manlngful,  an  examination  of  the 
nondlmenalonal  forma  nay  be  uaeful. 


According  to  reference  (n),  X^  ■ 2w“  " CDff)  • For  flight  condition  #1, 

■ 8.469  X 10'4.  Some  uncertainty  axlata  In  the  determination  of  trim 
due  to  prevlouely  dlacuaaed  dlfflcultlea  with  angle  of  attack  reference,  but 


(C  ) - la  probably  aufflclent  for  the  purposes  of  this  portion  of 

trl“  pU^S 

the  Investigation.  Table  XIV  lists  X^  reduced  to  nondlmenslonal  form,  and 
CD  extracted  from  the  combined  form,  for  all  "final11  runs  of  short  time 


length.  Unfortunately,  It  seems  clear  that  neither  (C^  - Cq^)  nor  Cq^  are 
much  closer  to  being  consistent  than  X Itself.  Nor  la  Cq  a consistent 


function  of 


trim* 


Reference  (t)  found  that  a of  0.6  corresponds  to  an  <0,)trl.  of 

approximately  6°,  and  that  a Cq  of  0.287  was  Indicated  for  this  condition 

Of 

(without  consideration  of  drag  due  to  elevator  required  for  trim).  Thla  value 
of  Cq  la  considerably  smaller  than  those  determined  by  the  identification 


program.  Reference  (t)  does  Indicate  that  Cq  Increases  rapidly  with  Increasing 


(o)  trim*  h*vln«  * value  of  °*573  •*  (^trlm  " 10  » and  1*123  #t  (or)trljI 


12. 5#. 


INTERPRETATION  OF  Mj  RESULTS 

* 2 

The  equation  (reference  (n))  M«  ■ C*  can  be  used  to  find  the 

» •iy 

Cja  equivalents  of  the  Identified  ! $ values.  Table  XV  lists  the  results 
^e  * 

of  these  calculations. 
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TA3LE  XIV 

FLIGHT  CONDITION  RESULTS  FOR  CDfy 


Run 

X 

a 

<CL  " %> 

<CL>trim 

Cn 

A Priori 

-23.23 

-0.492 

0.683 

1.175 

DR3-1W2 

-5.7 

-0.111 

0.625 

0.736 

DR3-1W1DEZ 

+0.4 

+0.008 

0.625 

0.617 

DR2-IW3 

-10.92 

-0.182 

0.535 

0.717 

DR4-1W4 

-10.63 

-0.196 

0.593 

0.789 

DR5-1W3 

+5.11 

+0.092 

0.582 

0.490 

DR6-1W3S 

-8.55 

-0.159 

0.601 

0.760 
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TABLE  XV 

FLIGHT  CONDITION  h RESULTS  FOR  Cm*e 


Run 

M‘e 

A Priori 

-9.68 

-1.18 

DR3-1W2 

-7.20 

-0.803 

DR2-1U3 

-9.32 

-0.891 

DR4-1W4 

-8.26 

-0.875 

DR5-1W3 

-8.65 

-0.899 

DR6-1W3S 

-8.15 

-0.875 
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Reduction  to  nondimen3ional  forms  tends  to  reduce  the  differences  between 
the  identified  parameters.  The  smallest  Cm)t  is  only  an  11%  reduction  from 

e 

the  largest  identified  Cm>e,  while  the  smallest  M*e  was  a 23%  decrease  from 

the  largest  identified  M«  . However,  observation  of  the  Cm#  results  reveals 

e e 

that  the  identified  results  are  all  quite  different  from  the  a priori  value. 

The  Cm*  determined  by  the  wind  tunnel  tests  reported  in  reference  (t)  is 

e 

1.29  for  *e  between  -10°  and  -15*,  which  is  slightly  higher  than  the  a priori 
value  and  much  different  from  the  identified  parameters.  However,  the  wind 
tunnel  tests  also  indicate  that  elevator  effectiveness  drops  off  rapidly  with 
increasingly  negative  *e  beyond  *e  - -15*,  suggesting  that  a linear 

e 

might  not  be  appropriate.  Yet  the  reference  (t)  wind  tunnel  data  are  for  a 
no-flaps  configuration,  while  the  flight  data  were  taken  for  a half-flaps 
configuration.  The  deflection  of  the  flaps  certainly  could  change  the  elevator 
effectiveness  and  the  elevator  angle  at  which  effectiveness  decreases.  Thus, 
the  comparison  of  the  reference  (t)  data  and  the  flight  results  may  not  be 
meaningful,  but  the  a priori  information  was  generated  for  the  half-fla|  case 
and  no  simple  explanation  exists  for  the  discrepancy  between  this  value  nd 
che  flight  results. 

SUMKVRY  OF  FLIGHT  CONDITION  #2 

Identification  results  for  flight  condition  #2,  reduced  with  respect  to  an 
assumed  known  M-,  are  given  in  Table  771.  In  general,  the  results  are  consistent 
within  themselves , although  not  necessarily  in  agreement  with  the  a prio-i 
values.  Unreduced  results  are  given  in  Appendix  C. 

DMA  PUN  19 

Data  run  19  featured  a pseudo-white  noise  elevator  input  with  appro;: mutely 
zero  mean.  The  continuous  reversal  movements  of  the  control  allowed  th  . onlroi 
deflection  magnitude  to  be  relatively  large  without  causing  any  motion 
variables  to  exceed  the  range  for  which  linear  dynamics  may  reasonably  be 
considered.  Thus  the  signal  to  noise  ratio  for  the  elevator  signal  was  high 
enough  to  permit  processing  without  first  smoothing  the  control  data.  The 
resultant  time  history  matches  are  shown  in  Figure  17.  General  agreement: 
seems  evident  although  the  matches  are  by  no  means  perrect.  The  largest 
discrepancy  occurs  due  to  an  elevator  noise  spike  at  approximately  four  seconds 
The  angle  of  attack  match  has  a persistent  offset  problem,  however,  and  the 
normal  acceleration  response  of  the  model  seems  to  consistently  lag  the 
flight  data.  The  pitch  acceleration  response  of  the  model  is  too  noisy  to 
permit  any  conclusions. 

DMA  RUN  10 

The  time  history  matches  for  data  run  10  are  shown  in  Figure  18.  The 
signa L/noise  ratio  for  the  elevator  data  (solid  line)  is  clearly  poor.  Shown 
dotted  on  the  elevator  trace  is  the  smoothed  elevator  time  history  actually 
used  as  input  r.o  the  identification  program.  The  time  history  match  is 
lather  good,  bit  the  maneuver  is  so  small  that  it  would  seem  dangerous  to 
place  too  much  weight  on  the  results. 
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TABLE  XVI 

FLIGHT  CONDITION  #2  FINAL  RESULTS 
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DATA  RUN  14 

Only  IS  seconds  of  data  were  available  in  data  mn  14.  The  tiae  history 
matches  obtained  for  the  shortened  run  are  shown  In  Figure  19.  The  magnitudes 
of  the  peaks  of  the  angle  of  attack  and  pitch  angle  responses  In  their  positive 
directions  are  significantly  larger  for  the  model  than  for  the  flight  data. 

Alec,  the  model  normal  acceleration  response  again  seems  to  lag  the  flight  data. 

Although  the  results  for  Individual  stability  derivatives  will  be  discussed 
in  later  sections,  It  Is  worth  noting  here  that  the  largest  values  of  Z , M , 

^ (V 

and  and  the  smallest  value  of  Zj#,  were  obtained  from  this  data  run.  In 

particular,  the  Z value  of  -1656.  is  the  only  Identified  result  approaching 
o 

the  a priori  value  of  -1669. 

DATA  RUN  13 

The  control  Input  shape,  and  the  errors  in  the  resulting  time  history 
natch,  for  data  run  13  are  nearly  identical  to  data  run  14.  Yet  the 
Identified  Z for  data  run  13  Is  191  lower  than  that  obtained  from  data  run  14, 

ft 

and  X , M , and  M are  also  lower.  The  time  hintory  comparisons  for  data  run 

rt  n q 

13  are  shown  In  Figure  20. 

DATA  RUN  18 

Tims  history  comparisons  for  data  run  18  are  shown  in  Figure  21.  The 
modal  responses  for  angle  of  attack  and  pitch  angle  roughly  correct  in 
shape  but  have  a substantial  bias  In  the  direction  of  positive  response.  This 
phenomenon  is  not  a simple  shift  of  the  entire  reuponse  curves,  because  the 
model  response  curves  and  the  flight  data  are  similar  at  both  the  beginning 
and  end  of  the  data  run.  The  normal  acceleration  shows  that  the  model  response 
lags  the  flight  data  "id  the  amplitude  of  the  model  response  Is  Insufficient 
in  the  positive  dlrec  ion  to  match  the  flight  data. 

The  ax  time  history  comparison  In  this  data  run,  as  veil  as  all  the 

previously  discussed  data  runs  for  flight  condition  #2,  contains  little 
information  because  of  the  high  noise  level  In  the  a^  flight  data  and  the 

relatively  low  level  of  excitation.  It  Is  quite  significant,  however,  that 
the  model  and  flight  responses  are  In  the  same  direction.  It  should  be 
recalled  that  the  responses  tended  to  be  In  opposite  directions  in  flight 
condition  #1  time  history  comparisons,  leading  to  speculation  that  a sign 
error  erlsted  In  either  the  model  or  the  measurement  system.  The  somewhat 
larger  ax  excitations  In  the  data  runs  of  flight  condition  #2  reveal  that  the 

signs  are  correct  in  general,  although  the  a^  matching  problems  in  flight 

condition  #1  'nay  indicate  s deficiency  in  the  model  for  that  configuration. 

The  results  from  data  run  18  feature  the  sswllest  7.  , M , and  M , and 

v ft  q 

the  largest  value  of  Zg#,  of  any  flight  condition  #2  data  run.  This  situation 
is  exactly  opposite  that  of  data  tun  14. 
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FIGURE  19.  TIME  HISTORY  COMPARISON  DR14-1W3S2SH 
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FIGURE  ZD.  THE  HISTORY  CGME&RISON  DR13-1W2S2 


NA  DC-74181-30 


i ' 


83 


FIGURE  20  (CONT).  TIME  HISTORY  COMPARISON  DR13-1W2S2 


FIGURE  21.  TIME  HISTORY  COMPARISON  DR18-1W3S2 
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FIGURE  21  (CONT).  TIME  HISTORY  COMPARISON  DR18-1W3S2 
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DATA  RUN  11 

Figure  22  is  a time  history  comparison  for  data  run  11.  Note  that  even 
with  two  iterations  of  smoothing,  the  elevator  noise  spikes  cause  considerable 
model  response  which  is  not  present  in  the  flight  data.  Figure  23  is  a similar 
comparison  for  data  run  11  with  elevator  set  to  zero  following  the  significant 
input.  Once  again  the  angle  of  attack,  pitch  angle,  and  normal  acceleration 
model  responses  overshoot  the  flight  data  in  the  time  inmediateiy  following 
the  input.  The  identification  attempts  not  including  elevator  zeroing  were  not 
considered  "final",  and  the  results  were  not  used.  The  time  history 
comparison  plots  were  shown  to  illustrate  the  severity  of  the  elevator  noise 
problem. 

Unfortunately,  a plot  of  che  70  second  identification  processing  of  data 
run  11  was  not  available  due  to  computer  malfunction.  Airspeed,  the  quantity 
of  particular  interest,  varied  to  -25  ft/sec  and  returned  to  approximately 
the  original  trim  value  during  the  data  run.  The  stability  derivative  results 
(with  all  but  the  speed  derivatives  held  at  the  values  obtained  in  DR11-1W2S2DEZ) 
are  shown  in  Table  XVI. 

DATA  RUN  12 

The  pilot  input  for  this  data  run  was  insufficient  to  cause  motions  large 
enough  to  permit  identification  of  the  stability  derivatives  principally 
involved  in  the  short  period  motion.  The  intent  of  the  input  was  the  generation 
of  phugoid  motions,  in  which  it  succeeded.  An  87.5  second  record  length 
segment  of  data  run  12  contained  11;  cycles  of  phugoid  motions,  including 
velocity  variations  between  +41  ft/sec  and  -30  ft/sec. 

Once  again,  a plot  of  the  time  history  comparison  was  unavailable,  but 
the  stability  derivative  results  are  shown  in  Table  XVI.  Start-up  values  for 
DR12-U5W3S2DEZ  were  not  the  usual  a priori  values,  but  the  final  results  of 
DRU-U4W2S2DEZ. 

STABILITY  DERIVATIVE  RESULTS 

The  variation  in  identified  stability  derivatives  from  run-to-run,  and 

with  respect  to  their  a priori  values,  was  similar  to  flight  condition  #1 

with  certain  exceptions.  As  shown  in  Table  XVI,  results  for  X were  somewhat 

o 

more  consistent,  but  still  covered  a wide  range.  Both  (C. ) . and  Cn 

L trim 

would  be  expected  to  be  quite  small  for  flight  condition  #2,  and  any  measure- 
ment or  the  differences  between  them  could  certainly  be  subject  to  considerable 
error.  The  a priori  value,  which  is  larger  in  magnitude  than  any  identified 
and  opposite  in  sign  to  all  of  the  identified  X^  results,  could  also  be 

subject  to  error  for  the  same  reason.  Calculation  of  the  a priori  value 
required  the  measurement  of  the  slope  of  a very  shallow  vs.  * curve. 

The  results  for  M$e  were  remarkably  consistent.  This  phenomenon!  may 
indicate  that  the  identification  problems  were  unique  to  flight  condition 
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FIGURE  22.  TIME  HISTORY  COMPARISON  DR11-IW2S2 
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#1,  or  just  the  opposite  could  be  true  if  the  algorithm  were  in  fact  simply 
retaining  the  a priori  value  due  to  a lack  of  sufficient  information. 

The  Mq  identification  results  are  rather  consistent,  but  quite  different 

from  the  a priori  value,  the  latter  only  having  about  half  the  magnitude  of 
the  former. 


CONCLUSIONS 


GENERAL 

1.  A modified  Nevton-Raphson  parameter  identification  technique  was  acquired. 

This  technique  was  shown  to  extract  nearly  perfect  sets  of  stability 
derivatives  from  computer-generated  data. 

2.  Carefully  documented  sets  of  flight  data  have  been  collected  and  provided 

to  several  contractors  and  agencies  for  the  purpose  of  parameter  identification. 
These  data  are  especially  suitable  for  parameter  identification  in  the  sense 
of  the  completeness  of  the  measurements  and  the  variety  of  selected  pilot 
inputs.  However,  the  signal-to-noise  ratios  for  many  of  the  measured 
quantities  are  not  high  enough  for  the  data  to  be  suitable  as  a test  case  for 
unproven  identification  techniques  or  for  automatic  identification  without 
human  involvement.  The  most  serious  problem  is  the  high  noise  level  on  the 
control  input  measurement. 

3.  The  measurement  and  recording  systems  currently  available  at  NAVAIRDEVCEN 
for  the  acquisition  of  flight  data  for  parameter  identification  are  barely 
adequate  for  the  purpose. 

4.  A great  deal  of  valuable  experience  has  been  gained,  and  expertise 
developed  at  NAVAIRDEVCEN,  in  the  acquisition  of  flight  data  for  parameter 
identification  purposes  and  in  the  use  of  parameter  identification  techniques. 

5.  No  significant,  consistent  trends  in  results  were  observed  as  a function 
of  pilot  input  type.  Without  any  type  of  quantitative  indication  of 
identification  accuracy,  it  is  not  possible  to  reach  any  conclusions  from  the 
current  analysis  concerning  the  effect  of  input  type  on  identification  performance. 

6.  Measurement  of  pitch  angular  acceleration  directly  with  a Statham  angular 
accelometer  is  a far  superior  method  to  the  calculation  of  pitch  acceleration 
based  on  nose  and  tail  normal  acceleration  measurements. 

FLIGHT  CONDITION  #1 

1.  Due  to  problems  with  the  performance  ot  the  airspeed  transducers,  filtering 
of  the  airspeed  measurement  was  necessary. 

2.  Proper  use  of  the  weighting  matrix  in  the  identification  program,  causing 
the  weighted  contributions  to  the  fit  error  by  the  various  time  history  matches 
to  be  equal,  was  essential  to  the  performance  of  the  identification  technique. 
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3.  Setting  the  value  of  the  control  input  to  zero  following  the  significant 
pilot  control  input  has  a significant  effect  on  the  identified  results,  but 

it  is  not  certain  that  this  technique  produces  valid  results,  since  information 
is  destroyed  as  well  as  noise. 

4.  The  algorithm  would  not  converge  unless  the  a priori  stability  derivatives 
were  reasonable  estimates.  In  the  absence  of  sufficient  information  in  the 
time  histories,  the  identified  results  are  heavily  dependent  upon  the  a priori 
values . 


5.  Only  one  data  run  (#3)  was  sufficiently  long  to  permit  a reasonable 
attempt  at  Identification  of  the  speed  derivatives.  Unfortunately,  these 
data  included  significant  atmospheric  turbulence,  which  was  not  included  in 
the  model  and  therefore  must  have  had  an  adverse  effect  on  the  identification 
process. 

6.  Based  on  identification  results,  the  best  estimate  of  X is  -0.046,  which 
agrees  with  the  a priori  value. 


7.  The  stability  derivative  was  not  consistently  identifiable.  Reduction 

of  results  to  Cn  still  does  not  reveal  any  single  value  or  function  of  (C  ) 

ft  L trim 

All  identified  values  were  substantially  smaller  than  the  a priori  value. 

8.  The  best  estimate  of  Z^  is  -0.39,  while  the  a priori  value  was  -0.27. 

9.  Based  on  best  results  from  data  run  segments  of  17.5  seconds  or  less,  the 
mean  estimate  for  Z was  -240.3  (compared  to  an  a priori  value  of  -233.8)  and 

ft 

the  standard  deviation  among  the  six  numbers  was  38.1. 

10.  The  best  estimate  of  My  is  -0.004,  compared  to  the  a priori  value  of  zero. 

11.  The  mean  estimate,  based  on  appropriate  results,  for  was  -5.11, 
exactly  identical  to  the  a priori  value,  with  a standard  deviation  of  0.34. 

12.  For  M , the  mean  estimate  was  -1.32  and  the  standard  deviation  was  0.1. 

q 

The  a priori  value  of  M was  -1.42. 

q 

13.  The  mean  estimate  of  Zg  was  -22.8  and  the  standard  deviation  was  1.81. 

The  a priori  value  was  -24. 

14.  The  mean  estimate  of  M*e  was  -8.22,  with  a standard  deviation  of  0.74. 

The  a priori  estimate  was  -9.68.  Although  there  was  less  scatter  in  the  Mge 
estimates  than,  for  Instance,  Z estimates,  the  large  discrepancy  between  the 

ft 

mean  estimate  and  the  a priori  value  was  a source  of  concern.  Reduction  of 

the  results  to  the  nondimensional  form  C,™.  , however,  led  to  an  increase  in 

m c 

the  consistency  of  the  identified  results  and  an  increase  in  the  difference 
between  the  mean  estimate  and  the  a priori  value. 
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15.  The  time  history  matches  were  generally  good.  The  most  serious  matching 
problems  occurred  in  the  angle  of  attack  time  histories. 

FLIGHT  CONDITION  #2 

1.  The  conclusions  reached  for  flight  condition  #1  concerning  airspeed 
signal  filtering,  use  of  weighting  matrix,  and  zeroing  of  control  input 
signal  are  also  applicable  to  flight  condition  #2. 

2.  The  signal  to  noise  ratio  in  the  elevator  signal  was  so  low  that  a smoothing 
process  was  required  before  the  elevator  signal  could  reasonably  be  used  in  the 
identification  program. 

3.  A priori  value  of  was  -0.016.  Only  two  estimates  were  obtained, 

-0.038  and  -0.099,  and  insufficient  agreement  exists  for  selection  of  a 
best  estimate. 

4.  A priori  value  of  was  -O.I22.  Estimates  were  -0.011  and  -0.233,  again 

insufficient  for  selection  of  a reasonable  estimate. 

5.  Estimates  of  M were  0.0016  and  0.0018,  which  agree  very  well  with  an  a 
priori  value  of  0.0017. 

6.  The  mean  estimate  of  was  20.5,  with  a standard  deviation  of  9.97 

(497.  of  the  mean).  Indicating  that  little  confidence  should  be  placed  in  this 
estimate.  The  a priori  value  of  -36.9  is  in  strong  disagreement  with  the 
identified  results. 

7.  The  mean  estimate  of  Zw  was  -1371.6,  in  contrast  with  an  a priori  value 

of  -1669.2.  This  discrepancy  seems  large  for  a derivative  which  should  be 
easily  identifiable.  The  standard  deviation  among  the  six  appropriate 
estimates  of  Z was  170.4  (127.  of  the  mean). 

8.  The  mean  estimate  of  ML  was  -22.4,  which  agrees  very  well  with  the  a 

Or 

priori  value  of  -22.5.  The  standard  deviation  was  less  than  27.  of  the  mean, 
indicating  very  consistent  identification. 

9.  The  mean  estimate  of  was  -7.33,  and  the  standard  deviation  of  the  run- 

to-run  values  was  0.67,  or  97.  of  the  mean,  indicating  fairly  consistent  results. 
Yet  the  a priori  value  was  only  -3.39  and  would  seem  to  be  seriously  in  error. 

10.  The  mean  estimate  of  Z*  was  -134.5,  in  contrast  to  the  a priori  value  of 

e 

-121.3.  The  standard  deviation  was  6.8  (57.  of  the  mean). 

11.  The  mean  estimate  of  M*  was  -49.3  and  the  standard  deviation  was  only 

e 

0.6,  or  about  17..  This  remarkably  consistent  set  of  identified  results  is  in 
reasonable  agreement  with  the  a priori  value  of  -52.1. 
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1?..  The  time  history  matches  were  only  fair.  A particular  pattern  of  excessive 
amplitude  in  the  model  response  for  angle  of  attack,  pitch  anple,  and  normal 
Acceleration  was  obtained  several  times. 
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APPENDIX  A 

LISTING  OF  PROGRAM  MASSAGE 
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M + M. 

q ~ 


U + Z 
o 

U -Z 


O r* 
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M*e  + 


Initial  Conditions:  X , Z , M 

O O O 


For  run  number  coding  refer  to  page  52. 
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Run  No. 


Derivative 

Name 

Start- 

Up 

DR2-1WL 

DR2-1W2 

DR2-1W3 

X 

u 

-.046 

-.066 

-.074 

-.076 

X' 

-.087 

-.049 

-.046 

-.041 

z 

u 

-.27 

.040 

.062 

.069 

Z' 

n 

-.877 

-1.142 

-1.16 

-1.17 

M' 

u 

1.59 

.28 

.23 

.21 

H' 

rv 

-4.59 

-4.95 

-4.99 

-5.02 

M' 

q 

-1.95 

-1.74 

-1.76 

-1.77 

-.092 

-.090 

N 

00 

O 

• 

t 

-.081 

e 

-9.63 

-9.16 

-9.24 

-9.28 

e 

X 

o 

- 

-.2  X 10“2 

-.3  X 10"2 

-.3  X 10' 

Z 

o 

- 

.8  X 10-2 

1.  X io*2 

.1  X 10' 

M 

o 

- 

-.2  X 10"1 

-.2  X 10-1 

-.2  X 10' 

Run  Time 

- 

17.5 

17.5 

17.5 

Trim 

Velocity 

- 

266.5 

266.5 

266.5 

C-4 


NA DC-74181-30 


TABLE  C-I  (CONTINUED) 

IDENTIFICATION  RESULTS  ^FLICHTCONDITION  ill) 


Run  No. 


Derivative 

Name 

Start- 

_JiE_ 

DR3-1 

DR3-1W1 

DR3-1W2 

X 

u 

-.046 

-.046 

-.028 

-.046 

X' 

rv 

-.087 

-.087 

-.047 

-.023 

z 

u 

-.27 

-.27 

-.26 

-.28 

Z1 

<v 

1 

• 

00 

-.877 

1 

• 

o 

w 

-.954 

M* 

u 

1.59 

1.60 

.935 

-.020 

M' 

-4.59 

-4.59 

-4.54 

-4.46 

M' 

q 

-1.95 

-1.95 

-1.97 

-1.75 

e 

-.092 

-.092 

-.089 

-.088 

M»« 

-9.63 

-9.61 

-7.56 

-7.15 

X 

o 

- 

-0.8  X 10-3 

-0.3  X 10"2 

-0.2  X 10 

Z 

o 

- 

0.4  X 10"3 

0.1  X io_1 

0.2  X 10 

M 

o 

- 

-0.6  X 10“4 

-0.2  X lO”1 

-0.2  X 10 

Run  Time 

- 

17.5 

17.5 

17.5 

Trim 

Velocity 

- 

246.7 

246.7 

246.7 

C-5 


NA DC-74181-30 


TA3LE  C-I  (CONTINUED) 
IDENTIFICATION  RESULTS  (FLIGHT  CONDITION 


Run  No, 


Derivative 

Name 

Start- 

Up 

DR3-2W 

DR3-2W2 

DR3-U2W1 

X 

u 

-.046 

-.046 

-.047 

-.054 

X' 

n 

-.087 

-.048 

-.030 

-.023 

Z 

u 

-.27 

-.356 

-.373 

-.414 

Z' 
r * 

-.877 

-.831 

-.817 

-.954 

M' 

u 

1.59 

-.57 

-.61 

-.73 

M' 

n 

-4.59 

-4.45 

-4.34 

-4.46 

-1.95 

-2.13 

-2.19 

-1.75 

Z\ 

-.092 

-.090 

-.089 

-.090 

e 

. 

"5. 

-9.63 

-7.07 

-7.10 

-7.15 

X 

o 

- 

-0.2  X 10’2 

-0.3  X 10-2 

-0.2  X 10' 

Z 

o 

mm 

0.2  X 10’1 

0.2  X 10"1 

0.2  X 10 

M 

o 

- 

-0.1  X io~l 

-0.2  X iO*1 

-0.2  X 10" 

Run  Time 

- 

35 

35 

35 

Trim 

Velocity 

- 

246.7 

246.7 

246.7 

NA  DC -74181-30 


TABLE  C-I  (CONTINUED) 

IDENTIFICATION  RESULTS  (FLIGHT  CONDITION  #1] 


Run  No. 


Derivative 

Name 

Start- 

Up 

DR4-1W1 

DR4-1W2 

DR4-1W3 

DR4-1W4 

X 

u 

-.046 

-.023 

-.014 

-.020 

-.036 

X' 

rv 

-.087 

-.057 

-.038 

-.043 

-.042 

z 

u 

CM 

• 

• 

-.16 

0 

CM 

• 

1 

CM 

CM 

• 

I 

CM 

CM 

• 

1 

Z' 

(V 

-.877 

-.922 

-.904 

-.887 

-.879 

M' 

u 

1.59 

-.67 

-.83 

-.93 

-.96 

M* 
r » 

-4.59 

-4.30 

-4.27 

-4.24 

-4.22 

S 

-1.95 

-1.72 

-1.76 

-1.79 

-1.81 

z; 

-.092 

-.096 

-.088 

-.090 

-.091 

e 

, 

\ 

-9.63 

-8.08 

-8.11 

-8.18 

-8.21 

X 

o 

• 

.2  X 10”2 

.9  X 10"3 

-.1  x io"3 

-.7  X 10 

Z 

o 

- 

-.4  X 10'3 

-3 

-l.  X 10  J 

-.2  X 10-2 

-.2  X 10' 

M 

o 

- 

.6  X 10‘2 

.3  X 10-2 

-.4  X 10“3 

-.3  X 10’ 

Run  Time 

17.5 

17.5 

17.5 

17.5 

Trio 

Velocity 


253.2  253.2  253.2  253.2 


NA DC-74181-30 


TABLE  ^-1  (CONTINUED) 

IDENTIFICATION  RE~SULTS  (fLICHTCONDITION  #1) 


Run  No. 


C-9 


Derivative 

Name 

Start- 

Up 

DR5-1W1 

DR5-1W2 

DR5-1W3 

X 

u 

-.046 

.034 

.035 

.034 

X' 

rv 

-.087 

.023 

.024 

.020 

z 

u 

j 

-.27 

-.35 

-.36 

-.35 

Z' 

(V 

-.877 

-.919 

-.922 

-.941 

M' 

u 

1.59 

1.29 

1.29 

1.25 

M1 

r* 

-4.59 

-4.80 

-4.81 

-4.86 

-1.95 

-1.87 

-1.88 

-1.85 

z\ 

-.092 

-.113 

-.113 

1 

• 

o 

W 

e 

m; 

-9.63 

-8.66 

-8.67 

-8.60 

e 

X 

o 

.7  X 10"3 

.7  X 10-3 

.7  X IO-3 

z 

o 

- 

.60  X 10"2 

-2 

.6  X 10 

.6  X 10‘2 

M 

o 

- 

-1.  X io’2 

-1.  x io"2 

-.1  x io-1 

Run  Time 

- 

17.5 

17.5 

17.5 

Trim 

Velocity 

- 

255.7 

255.7 

255.7 

NA DC-74181-30 


TABLE  C-I  (CONTINUED) 

IDENTIFICATION  RESULTS  (FLIGHT  CONDITION  #1) 


Run  No 

• 

Derivative 

Name 

Start- 

Up 

DR6-1W1 

DR6-1W2 

DR6-1W2SH 

DR6-1W3SH 

X 

u 

-.046 

.024 

.025 

.018 

0 

X' 

n 

-.087 

-.052 

-.055 

-.044 

-.034 

Z 

u 

-.27 

-.38 

-.38 

-.25 

-.24 

Z‘ 

n 

-.877 

-.941 

-.945 

-.937 

-.927 

M' 

u 

1.59 

.56 

.54 

.59 

.47 

M' 

rv 

-4.59 

-4.51 

-4.53 

-4.54 

-4.41 

-1.95 

-1.76 

-1.79 

-1.81 

-1.86 

• 

-.092 

-.094 

• 

• 

o 

VO 

o\ 

-.089 

• 

• 

o 

00 

Ov 

«; 

-9.63 

-6.74 

-6.81 

-7.97 

-8.10 

e 

X 

o 

• 

.2  X 10‘2 

.2  X 10”2 

.2  X 10”2 

.2  X 10‘2 

Z 

o 

m 

-.1  X 10*1 

-.1  X 10’1 

-.1  X io"1 

-.1  X io"1 

M 

o 

- 

.1  X 10-1 

.2  X 10"1 

.2  X 10"1 

.2  X 10“l 

Run  Time 

17.5 

17.5 

10.0 

10.0 

Trim 

Velocity 

251.6 

251.6 

251.6 

251.6 

NA  DC-74181-30 


1A3LE  C-II 


IDENTIFICATION 

RESULTS  (FLIGHT  CONDITION  #2) 

Run  No. 

Derivative 

Name 

Start- 

Up 

DR10-1W1S1 

DR10-1W2S1 

DR10-1W2S2 

X 

u 

-.016 

-.0034 

.0017 

.0010 

X* 

iy 

-.054 

.0065 

.011 

.010 

z 

u 

-.12 

+.179 

+.156 

+.155 

Z’ 

tv 

-2.46 

-1.89 

-2.05 

-2.06 

M’ 

u 

1.34 

-.375 

-.63 

-.63 

M' 

<v 

-18.6 

-19.07 

-19.0 

-19.0 

-4.96 

-8.59 

-8.15 

-8.04 

zl 

00 

H 

• 

1 

-.19 

-.19 

-.19 

e 

K 

-51.8 

-49.0 

-49.4 

-49.8 

e 

X 

o 

«a 

.15  X 10“3 

.3.X  10"3 

.3  X 10”3 

Z 

o 

-2 

-.5  X 10 

-.5  X 10"2 

-.5  X 10'2 

M 

o 

m 

.6  X 10"3 

.5  X 10"2 

-2 

.5  X 10 

Run  Time 

- 

17.5 

17.5 

17.5 

Trim 

Velocity 

- 

683.7 

683.7 

683.7 

i 


C-ll 


NA DC-74181-30 


TABLE  C-II  (CONTINUED) 
IDENTIFICATION  RESULTS  (FLIGHT "CONDITION  #2) 


Run  No. 


Derivative 

Name 

Start- 

DR11-1W1S2 

DR11-1W2S2 

DR11-1W1S2DEZ 

X 

-.016 

-.026 

-.038 

i 

• 

o 

o 

u 

X' 

-.054 

.012 

.028 

.020 

rx 

Z 

u 

-.12 

.029 

.092 

.172 

Z' 

-2.46 

-2.28 

-2.13 

-1.80 

a 

M' 

1.34 

CO 

a 

.751 

.32 

u 

M' 

-18.6 

-19.0 

-19.1 

-19.1 

% 

-4.96 

-9.78 

-9.05 

-7.29 

K 

e 

-.18 

00 

H 

• 

1 

-.19 

0 

CM 

a 

1 

K 

-51.8 

-48.8 

-49.1 

-50.3 

e 

X 

o 

m 

-.2  X 10-3 

-.5  X 10*3 

.2  X 10‘3 

Z 

o 

• 

.2  X 10"1 
. . -2 

.2  X 10"1 
. -2 

,1  x io"1 
-2 

M 

-.4  X 10 

-.7  X 10 

-.4  X 10 

o 

Run  Time 

17.5 

17.5 

17.5 

679.2  679.2 


Trim 

Velocity 


679.2 


fA  DC-74181-30 


TABLE  C-II  (CONTINUED) 
IDENTIFICATION  RESULTS  (FLIGHT  "CCTIDITIOH  »2) 


Run  No. 


Derivative  Start- 

Name  Up 


DR12-U5W1S2DEZ  DR12-U5W2S2DEZ  DR12-U5W3S2DEZ 


X 

u 

1 

• 

o 

O' 

-.038 

-.040 

-.038 

X' 

-.054 

-.054 

to 

o 

. 

1 

.017 

z 

u 

-.12 

1 

• 

N) 

O 

v4 

-.215 

-.233 

Z* 

f* 

-2.46 

-2.03 

-2.03 

-2.00 

M' 

u 

1.34 

1.40 

1.28 

1.57 

M' 

r* 

-18.6 

-22.0 

-22.0 

-19.3 

-4.96 

-7.0 

-7.0 

-8.69 

Zi 

-.18 

-.2 

-.2 

-.20 

e 

K 

-51.8 

-50. 

-50. 

-49.3 

e 

X 

o 

- 

-.2  X 10*2 

-.2  X 10"2 

-.1  X 10 

Z 

o 

- 

-.5  X 10"2 

-.5  X 10“2 

-.6  X 10 

M 

o 

- 

.8  X 10"2 

.1  X io*2 

.4  X 10 

Run  Time 

- 

87.5 

87.5 

87.5 

Trim 

Velocity 

- 

669.9 

669.9 

669.9 

fA  DC-74181-30 


TABLE  C-II  (CONTINUED) 

IDENTIF TrATTOM  RESULTS  fr LIGHT  CONDITION  #2) 


1 


Run  No. 


NA  DC-74181-30 


TABLE  C-1I  (CONTINUED) 
IDENTIFICATION  RESULTS  ^FLIGHT  CONDITION  »2) 


Run  No. 

Derivative 

Name 

Start- 

Up 

DR14-IW1S1SH 

DR14-1W2S2SH 

DR14-1W3S2SH 

X 

u 

-.016 

.044 

.052 

.053 

X' 

1 

• 

o 

i/I 

.043 

.0<+2 

.041 

Z 

u 

-.12 

+.111 

+.021 

1 

• 

o 

** 

Z' 

f* 

-2.46 

-2.02 

-2.28 

-2.41 

M' 

u 

1.34 

.40 

.73 

.91 

M' 

r» 

-18.6 

-19.5 

-19.3 

-19.2 

H' 

q 

-4.96 

-10.36 

-9.78 

-9.65 

K 

-.18 

0 

01 
• 
i 

-.19 

-.18 

e 

-51.8 

-47.7 

-48.6 

-48.7 

e 

X 

o 

- 

.2  X 10"2 

.2  X 10-2 

.2  X 10"2 

Z 

o 

- 

.1  X io"1 

.1  X io"1 

.1  x io-1 

M 

o 

- 

-.1  x io-1 

-.1  x io’1 

-.1  x io’1 

Run  Time 

* 

15.0 

15.0 

15.0 

Trim 

Velocity 

- 

687.2 

687.2 

687.2 

MDC-74181-30 


TABLE  C-II  (CONTINUED) 
IDENTIFICATION  RESULTS  (FLIGHT  CONDITION  #2) 

Run  No. 


Derivative 

Name 

Start- 

Up 

DR18-1W1S2 

DR18-1W2S2 

DR18-1W3S2 

X 

u 

H 

o 

• 

1 

+.044 

+.036 

.024 

X' 

rv 

-.054 

+.050 

.050 

.048 

z 

u 

-.12 

+.022 

-.021 

+.005 

Z' 

/V 

-2.46 

•1.66 

-1.44 

-1.64 

K 

1.34 

+1.94 

1.29 

1.00 

M' 

ft 

-18.6 

-21.0 

-19.7 

-19.3 

-4.96 

-9.23 

-8.88 

-8.62 

-.18 

-.21 

H 

CM 

• 

1 

-.21 

e 

-51.8 

-45.4 

-48.7 

-49.4 

e 

X 

o 

- 

.8  X 10"3 

.6  X 10’3 

.9  X 10"3 

Z 

o 

- 

.3  X 10"1 

.4  X 10"1 

.2  X 10"1 

M 

o 

m 

-.3  X 10’1 

-.9  X 10‘2 

-.5  X 10"2 

Run  Time 

m 

17.5 

17.5 

17.5 

Trim 

Velocity 

- 

682.8 

682.8 

682.8 

N*  DC-74181-30 


TABLE  C-II  (CONTINUED) 

IDENTIF ICATION  RESULTS  (FLIGHT  CONDITION  »2) 


Run  No. 


Derivative 

Name 

Start- 

Up 

DR19-1W 

DR19-1W2 

X 

u 

-.016 

-.097 

-.075 

X' 

rv 

-.054 

+.062 

+.036 

Z 

u 

-.12 

-.116 

-.013 

Z' 

n 

-2.46 

-2.00 

-1.98 

K 

1.34 

-.123 

+.276 

M' 

-18.6 

-20.0 

-19.55 

% 

-4.96 

-9.91 

-9.51 

zl 

-.18 

-.19 

0 

CM 

• 

1 

e 

t 

K 

-51.8 

-45.8 

-48.2 

e 

X 

o 

- 

-.2  X 10“2 

-.2  X 10" 

Z 

o 

- 

-.2  X 10"1 

-.3  X 10" 

M 

o 

m 

+.2  X iO"1 

+.8  X 10" 

Run  Time 

- 

17.5 

17.5 

Trim 

Velocity 

- 

690.0 

690.0 

Mote:  "S"  followed  by  digit  indicates  degree  of 

smoothing  for  elevator  signal 
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